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ABSTRACT 


For this thesis, an atmospheric propagation code named ANCHOR (Atmospheric NFS 
Code for High Energy Laser Optical pRopagation) was developed and utilized to study 
the propagation of high energy lasers in various atmospheric conditions and for numerous 
laser configurations. The ANCHOR code accesses existing industry databases to obtain 
relevant optical properties for various atmospheres and then uses scaling laws to simulate 
laser propagation through the defined environments. 

ANCHOR accounts for the effects of atmospheric diffraction, turbulence, 
platform jitter and thermal blooming on the laser beam, and outputs on-target irradiance 
and power-in-the-bucket profiles for a wide range of laser wavelengths. Several known 
physical trends associated with laser propagation will be reproduced, and the results will 
be compared to the industry accepted propagation code Wavetrain. 

The results of ANCHOR studies will indicate that the 100 kW-class high energy 
laser can effectively engage slow-moving targets at ranges greater than five kilometers in 
clear weather by delivering enough energy to melt 0.1 liters of one millimeter-thick 
aluminum aircraft skin in five seconds. For hazy, turbulent, and rainy conditions, the 
laser can effectively engage targets from ranges closer than three kilometers, but 
reasonable dwell times are only achieved for ranges closer than two kilometers. 
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I. 


INTRODUCTION 


The United States Navy is heavily invested in the researeh and development of 
high-energy lasers (HELs) onboard various surface and air platforms for offensive and 
defensive capabilities [1], [2], Such weapons, while extremely complex in design and 
expensive to manufacture, hold many advantages to currently-employed conventional 
weapons systems. From speed-of-light payload delivery, to an extremely deep magazine, 
to engagement flexibility for soft and hard kills, laser weapons are often touted as a 
significant part of the future of naval warfighting [3]. 

Unlike conventional weapons, however, laser weapons are uniquely affected by 
weather conditions of the environment in which they are employed [4], Effective 
development of these weapons is closely tied to the proper modeling and understanding 
of the role that atmospheric properties play in hindering or aiding laser propagation. 
Unfortunately, simply relying on experimental studies in order to fully parameterize laser 
behavior in various weather conditions is neither time-efficient nor cost-effective [4]. 
Physics-based computer modeling is therefore essential to guide HEE-based research. 

To this end, a streamlined and robust simulation code has been developed by the 
NFS Directed Energy Group. The code utilizes existing atmospheric databases and 
governing scaling equations to accurately and efficiently model laser propagation in 
various environments for numerous laser configurations. The code has been used to 
explore atmospheres with variable turbulence profiles, visibilities and precipitation, as 
well as laser effectiveness with various beam diameters, output powers, wavelengths, etc. 

The first part of this thesis will discuss the history, technology and advantages of 
directed energy weapons. Next, the thesis will explore the physical principles of 
atmospheric propagation, followed by atmospheric modeling, and then the simulation 
methods utilized to carry out the laser propagation studies. Easily, the thesis results will 
be discussed and analyzed, with conclusions tying together the quantitative and 
qualitative aspects of the research. 
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II. DIRECTED ENERGY WEAPON OVERVIEW 


A. DIRECTED ENERGY INTRODUCTION 

A directed energy (DE) weapon emits energy directed at a target in order to inflict 
damage without the use of a projectile [5]. The energy harnessed by a DE weapon is 
usually in the form of electromagnetic radiation, which includes microwave weapons and 
lasers, but also can include particle-beams and sonic weapons. Currently, DE weapons of 
military interest utilize a laser (light amplification by stimulated emission of radiation) to 
generate and deliver a concentrated spot of light at a specific wavelength and of high 
enough irradiance on the surface of a target to cause damage. The wavelength range of 
the laser light typically lies in the visible to infrared regime (~0.5 pm to 20 pm), and the 
energy propagates as a narrow beam at the speed of light to reach the target. Damage is 
primarily inflicted by the high-energy laser “by light dwelling on the target for an 
extended period of time to accumulate deposited energy and inflict thermal damage” [4]. 

Easer lethality, defined as “the capability of a weapon system to render a target 
non-functional,” [4] greatly depends on the target material the amount of power 
transmitted to the target. Thermal damage is considered to be the main lethality 
mechanism for lasers, resulting in the burning, melting, penetration or structural integrity 
failure of the target [4]. Eor typical targets, such as the carbon-carbon ceramic nosecone 
of a surface-to-air missile (SAM), an estimated 12.1 kilowatts (kW) of power over a 
dwell time of 4 seconds is required to ablate a depth of 1.0 cm on the target [6]. While a 
laser emitting such an output power level may not appear to be powerful enough to 
neutralize a missile traveling at Mach 3.0, a DE weapon of three-times the power and 
dwell time of the example may very well cause enough structural damage to “bring about 
catastrophic aerodynamic instability which would tear the missile apart” [6]. 

The first laser, constructed by Theodore H. Maiman in 1960, was a pulsed ruby 
laser with approximately one watt of output power [4]. Almost immediately following 
this initial invention, the U.S. Department of Defense (DOD) heavily invested into the 
research and development of lasers as operational military weapons systems. Several 
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military entities have since developed multiple HEL weapons systems for operational 
deployment. The U.S. Air Force has constructed the Advanced Tactical Laser and the 
Airborne Laser (ABL) test-bed, both of which employ aircraft-mounted systems utilizing 
chemical oxygen iodine lasers (COILs) with output powers of 100 kW and one megawatt 
(MW), respectively [4], 

The U.S. Navy has also developed several sophisticated weapons platforms, 
including the Mid-Infrared Advanced Chemical Laser (MIRACL), the Laser Weapon 
System (LAWS), and the Maritime Laser Demonstrator (MLD). The MIRACL system, 
first made operational in 1980, employs a deuterium-fluoride laser with ~one MW output 
power, while LAWS and MLD employ a solid-state fiber laser and a solid-state slab laser 
with output power of 32 kW and ~100 kW, respectively [4]. The latter laser system, the 
MLD, also became the first to be tested on and powered by a ship in 2011 onboard the 
U.S. Navy defense test ship USS Paul Foster (DD-964). The most significant 
achievement in the implementation of a DE weapon system onboard a naval platform, 
however, is the slated operational deployment of the aforementioned LAWS on the USS 
Ponce (LPD-15) in the Persian Gulf during summer 2014 [7]. 

B. ADVANTAGES OF HELS 

Conceptually, directed energy weapons systems have significant superiority to 
most conventional weapons systems. The main advantages of DE weapons lie in the 
precision and speed of payload delivery, low cost per engagement, and an essentially 
unlimited magazine [4]. A laser is capable of focusing to a spot several centimeters in 
diameter on a target from a range of up to 10 km [6]. Unlike most conventional weapons, 
such as explosive missiles and kinetic rounds, DE weapon precision also greatly reduces 
the possibility of collateral damage. Furthermore, due to the speed of light delivery of 
damage, a directed energy weapon has the potential of engaging and destroying a target 
soon after it is detected, whereas munitions from a conventional weapon must travel 
toward the target at a much slower rate and thus engage it at a close range. Radar may 
also be utilized for multi-target tracking, allowing for the continuous adjustment of the 
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HEL to maintain the beam on a single target while predicting the future location of any 
other targets for subsequent engagement. 

DE weapons systems are comprised of all-electric components, and complete 
integration into a naval ship’s power grid provides for a magazine potentially limited 
only by the ship’s power. All-electric DE weapons are thus by design a relatively cost- 
effective platform, and several second “shots” from an HEE are only theoretically tied to 
the price per gallon of the platform’s fuel. Although the initial overhead cost of DE 
weapon design and implementation has caused significant postponement in a budget- 
constrained environment, the cost per shot compared to an advanced tactical missile 
makes the HEE a worthy long-term investment. 

C. HEL TECHNOLOGY 

The primary types of DE weapons the U.S. Navy is considering for employment 
on expeditionary naval platforms are solid-state lasers (SSEs) and free-electron lasers 
(EEEs) [1]. Regardless of specific platform, the ideal HEE weapon system should exhibit 
several design traits. The system should have high wall-plug efficiency to generate the 
most optical power with as little input power as possible. The HEE should utilize a gain 
medium with high thermal conductivity for rapid heat removal from the system. The 
lasing medium should also be resilient against heat and/or dielectric breakdown to allow 
for scaling to higher power without permanent damage to the weapon. Einally, the output 
optical wavelength should reside in a region of atmospheric transparency, so as to 
minimize the amount of attenuation of the laser energy and to maximize the amount of 
energy delivered to the target. 

1. Solid-State Lasers 

SSL technology utilizes a solid substrate that holds the dopant in which the lasing 

transitions occur. Because this technology utilizes a material medium, substrate damage 

threshold limits output power [4]. SSLs have several advantages over EEEs, including 

their relatively compact and lightweight design (for low or modest power outputs), 

relative power efficiency, and relative operational safety in the absence of potential stray 

radiation. Most importantly, SSLs are more technically mature than EEEs, and can 
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potentially see more near term deployment availability. However, SSLs also have several 
disadvantages to FELs, ineluding the laek of sealability to MW-elass, a fixed optieal 
wavelength primarily determined by the gain medium, and diffieulty in obtaining 
aeeeptable beam quality for high-energy applieations [1], [4], Furthermore, DE 
applieations require the eombination of multiple solid-state lasers to attain the neeessary 
power output to infliet damage on target, whieh inereases their size and weight and 
degrades beam quality. 

The two types of SSEs whieh have been employed as taetieal weapons by the 
DOD are slab lasers and fiber lasers [1]. A slab laser eouples pump light with laser 
energy through one side of a reetangular medium slab, while removing heat from the 
opposite side of the slab. A fiber laser, on the other hand, utilizes a long, thin, eylindrieal 
rod medium, eovered by eladding, to eouple light to. While slab lasers have the 
geometrieal advantage of a larger aspeet ratio whieh allows for more light eoupling (as 
output power is proportional to aspeet ratio), fiber lasers are more lightweight, power 
effieient, durable and have eomparable output power to slab lasers with potentially better 
beam quality in a smaller deviee. 

Furthermore, optieal waveguides in fiber lasers span many meters in length, 
providing a mueh larger surfaee area that allows for easier eooling than a slab laser. 
While sueh advantages indieate that the fiber SSE is a theoretieally superior teehnology 
to the slab SSE, the latter has aetually demonstrated higher output power with the ~105 
kW of the MED, as opposed to the ~33 kW fiber SSE EAWS [1]. Ultimately, however, 
the most signifieant issue with SSEs, whieh the FEE may be able to address, lies in the 
diffieulty of eoherent beam shaping from eombining multiple SSEs, whieh is neeessary to 
attain higher output power [4]. 

2. Free-Electron Lasers 

Free-eleetron lasers utilize free (unbound) eleetrons to generate light [4]. An 
eleetron gun, eonsisting of an eleetron-generating eathode and an aeeelerator, ereates the 
initial eleetron beam. While several types of eathodes ean be employed to generate 
eleetron beams, thermionie, or “hof ’ eathodes are preferred for FEE applieations due to 
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their much higher power density [4], The electrons are accelerated to relativistic speeds 
within a linear accelerator, and subsequently pass through an undulator, consisting of a 
series of magnets which generate an alternating magnetic field, causing the electrons to 
“wiggle” and emit light. The interaction of the electrons with the light in the resonator 
produces stimulated, coherent emission. The wavelength of the light produced depends 
on the energy of the electrons and the properties of the undulator. Furthermore, wall-plug 
efficiency of the system can be increased by recovering most of the energy of spent 
electrons [4]. 

Because FELs utilize free electrons as the gain medium, elimination of waste heat 
is not an issue and thus their output power is thought to be scalable to megawatt class. 
Furthermore, the FEE optical wavelength is tunable, and thus can be chosen for optimal 
atmospheric transparency. An EEE can also obtain excellent beam quality at high output 
power for long-distance propagation and a relatively energy-efficient design. However, 
EEEs tend to be rather large in size (~ 20 meters long), heavy, and expensive to design 
[1]. EEEs also require shielding to contain radiation produced by stray electrons, and the 
overall technology is relatively less mature compared to SSEs. 

D. HEL MISSIONS AND PLATFORMS 

The missions of a DE weapon are largely dependent on its output power level. 
Directed energy weapons with lower power output (10-100 kW), such as the SSE, are 
better suited against slower moving targets, such as small watercraft and aircraft. Such 
DE weapons could immobilize boats and down aerial drones, as well as offer the 
potential of “soft kills” by blinding incoming missiles. Megawatt-class DE weapons, such 
as the EEE, will be primarily slated for cruise missile defense. The current EEE DOD 
project, the Innovative Naval Prototype, maintains a planned 100 kW of output power 
and is a stepping stone in the future transition to the MW-class devices. 

The HEE weapons systems of both kW and MW-class could be fittingly 
employed on several naval platforms, such as the future all-electric ships, littoral combat 
ships (ECS), guided-missile destroyers (DDG-51 and DDG-1000), and the future Eord 
class aircraft carriers, such as the Gerard R. Ford (CVN-78). Due to weight and size 
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constraints, only SSLs could be potentially employed on Naval aviation platforms sueh as 
F/A-18E/F Hornet, F-35 Joint Strike Fighter (JSF) [8], and Marine Corps helieopter 
platforms sueh as the UH-IY Huey and AH-IZ Cobra [6]. 
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III. ATMOSPHERIC PROPAGATION OVERVIEW 


The effectiveness of HEL weapons systems can be greatly affected by the 
medium through which the laser energy propagates. Because this medium is the Earth’s 
atmosphere, one needs to understand the various optical properties and physical 
principles associated with this environment to be able to develop systems that can operate 
effectively within it. 

A. ATMOSPHERIC EXTINCTION 

Molecules and aerosols suspended in the atmosphere absorb and scatter 
propagating laser energy as it travels to reach the target [4]. Absorption occurs when 
photons from incident radiation (such as laser light) transfer their energy to atmospheric 
constituents, thereby heating the local atmosphere. Scattering occurs when incident 
photons are redirected from the beam in random directions by atmospheric particulates. 

These combined effects are illustrated by Beer’s Eaw in Equation (1), which 
describes the transmission of light through a linear medium, following 

P(z) = P„e-’‘. ( 1 ) 

In Equation (1), Prepresents the total power delivered to a plane located a 
distance z from the source and is the initial output power at the beam director. The 

variable e represents the total extinction coefficient due to atmospheric absorption and 
scattering (in units of inverse length), and is further described by Equation (2) as 

^ = ( 2 ) 

In Equation (2), a refers to the absorption coefficients and refers to the 
scattering coefficients, while subscripts m and a denote the molecular and aerosol 
contributions, respectively. An increase in e results in greater attenuation rate of laser 
light by the atmosphere. 

I. Molecular Effects 

Table 1 displays a list of the major molecular constituents which make up the 
Earth’s atmosphere. 
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The Composition of the Atmosphere" 


Permanent constituents Variable constituents 


Constituent 

% by volume 

Constituent 

% by volume 

Nitrogen (N2) 

78.084 

Water vapor (H2O) 

0-0.04 

Oxygen (O2) 

20.948 

Ozone (O3) 

0-12 X 10-* 

Argon (Ar) 

0.934 

Sulfur dioxide (SO2)* 

0.001 X 10“-^ 

Carbon dioxide (CO2) 

0.036 

Nitrogen dioxide (NO2)* 

0.001 X 10--* 

Neon (Ne) 

18.18 X lO--* 

Ammonia (NH3)'' 

0.004 X 10-* 

Helium (He) 

5.24 X lO--* 

Nitric oxide (NO)* 

0.0005 X 10-* 

Krypton (Kr) 

1.14 X 10-^ 

Hydrogen sulfide (H2S)* 

0.00005 X 10- 

Xenon (Xe) 

0.089 X 10--^ 

Nitric acid vapor (HNO3) 

Trace 

Hydrogen (H2) 

0.5 X lO-* 

Chlorofluorocarbons 

Trace 

Methane (CH4) 

1.7 X 10-^ 

(CFCI3, CF2CI2 


Nitrous oxide (N2O)* 

0.3 X lO--’ 

CH3CCI3, ecu, etc.) 


Carbon monoxide (CO)* 

0.08 X lO--* 




“After the U.S. Standard Atmosphere (1976) with modifications. 
* Concentration near the earth’s surface. 


Table 1. Molecular composition of Earth’s atmosphere. From [9]. 


Ambient air is predominately composed of N 2 and O 2 , both of which are linear, 
homonuclear molecules with symmetric molecular bonds between each of their two 
atoms. Because these molecules only possess one normal vibrational mode that does not 
alter their overall dipole moments, propagating optical energy does not cause these 
molecules to vibrate. Thus, N 2 and O 2 also do not absorb electromagnetic radiation in the 
infrared spectral region [9]. 

Instead, other atmospheric constituents, such as H 2 O vapor and CO 2 , dominate the 
infrared atmospheric absorption spectrum. Even though these molecules contribute to a 
very small percentage of the total atmospheric composition, they are asymmetrically- 
shaped and exhibit multiple vibrational modes, making them strong absorbers of 
electromagnetic radiation [9]. The molecular absorption is extremely intricate and ideal 
laser performance requires that the spectral width of the beam fits within a wavelength 
window containing minimum absorption features. A typical molecular absorption 
spectrum for a wavelength range of 0.5 pm to 5.0 pm is shown in Figure 1. The gaps in 
the data at ranges ~2.5 to 3 pm, ~4.25 to 4.5 pm and ~4.75 to 5.0 pm are due to the 
absorption coefficient nearing infinity, which represents wavelength ranges of high 
absorption. 
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Molecular Absorption Spectrum 



Figure 1. Typical molecular absorption spectrum for the atmosphere. Gaps in 
the data correspond to high molecular absorption coefficient . 


Figure 2 displays a typical molecular scattering spectrum for the atmosphere over 
a wavelength range of 0.5 pm to 5.0 pm. Molecular scattering is described by an optical 
phenomenon known as Rayleigh scattering, in which shorter wavelengths of 
electromagnetic energy are scattered more readily than longer wavelengths. Rayleigh 
scattering will be further discussed in Chapter IV, Section B.4. 



Figure 2. Typical molecular scattering spectrum for the atmosphere. As a 
result of Rayleigh scattering (discussed in Chapter IV), shorter 
wavelengths of light scatter more readily than longer wavelengths. 
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Most of the photons scattered by molecules have the same wavelength as the 
incident photons; only their direction of propagation changes. This principle readily 
describes the phenomenon of the blue daytime sky and red sunsets. Because blue light 
more readily undergoes molecular scattering than red light, the daytime sky appears blue 
[9]. At sunset, the path length of light from the setting sun to the viewer through the 
atmosphere is longer, and by the time the light reaches the viewer’s eyes most of 
the shorter wavelength light has been scattered, leaving only the longer wavelength red 
light [9]. 

2. Aerosol Effects 

Absorption and scattering effects on light by aerosols are similar in nature to 
those by molecules. However, the extent of these effects depends greatly on aerosol size, 
shape and refractive index, the latter primarily defined by the material properties of the 
individual particles (water, dust, etc.). Figure 3 displays a plot of the aerosol absorption 
spectrum, while Figure 4 displays a plot of the aerosol scattering spectrum, both for a 
0.5 pm to 5.0 pm wavelength range. 

Aerosol Absorption Spectrum 

io-’r ^ 

lo-’i 

-I- ; ■ 

J 


10 ^ 



Figure 3. Typical aerosol absorption spectrum for the atmosphere. The peak in 
aerosol absorption at ~ 3 pm corresponds to presence of particulates 
which specifically absorb at those wavelengths. 
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Aerosol Scattering Spectrum 



2 2.5 3 3.5 

Wavelength (pm) 


X 10 


Figure 4. 


Typical aerosol scattering spectrum for the atmosphere. Aerosols 
scatter most wavelengths relatively equally. 


Because most aerosols are similar in size to the wavelength of infrared or visible 
light, they absorb most wavelengths approximately equally [5]. The notable increase in 
scattering at ~ 3 pm may be due to particulates such as sulfate (SO 2 ) or sea salt aerosols, 
which may absorb more at this specific wavelength because of their size. 

Similarly, aerosol particulates also scatter most wavelengths in the infrared and 
visible range relatively equally, which explains why clouds and haze appear 
approximately white to the human eye [9]. 

B. ATMOSPHERIC TURBULENCE 

Turbulence is a term generally describing a fluid flow regime within a medium in 
which the flow exhibits an irregular, chaotic and dissipative pattern, mostly in the form of 
randomly-varying eddies [10]. In the atmosphere, the turbulent flow of air is considered 
to be small-scale and primarily the result of convective air motion, vertical temperature 
differences, and wind shear. 

The major effects of turbulence on electromagnetic energy propagation stem from 
the spatial and temporal fluctuations in the refractive index of air [10]. The relationship 
between absolute temperature T (in K), wavelength X (in pm), pressure P (in mbar) 

and the refractive index of air n can be described by [10] 
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( 3 ) 


n ~ 1 + 77.6x10^® [l + 7.52xl0^'r"]-. 

Sufficient turbulence in the local atmosphere leads to variations in the local 
ambient temperature, which in turn leads to variations in the refractive index. For HEL 
weapons systems to function most effectively, the propagation of laser energy through 
the atmosphere requires that the light rays emitted by the laser maintain transverse 
coherence so that they constructively interfere to form a high irradiance spot size on 
target. However, refractive index variations due to turbulence cause different parts of the 
laser beam to experience different path lengths as the beam propagates, which disrupts 
constructive interference at the target and can even cause the beam to spread out and/or 
scintillate (i.e., breakup into beamlets) [10]. 


The Fried parameter (units of length) can be used to estimate the effects of 
turbulence on laser propagation. The parameter defines a circular diameter over which 
the laser beam maintains coherence (no significant change in phase) in the transverse 
direction throughout the propagation distance [4]. In other words, estimates the 
maximum transverse separation at the source between two individual light rays within 
which the two rays will constructively interfere at the target. A lower value implies 

stronger turbulence. If the value is significantly smaller than the beam director size, 
the laser beam will break up into many smaller, incoherently radiating beamlets. 


The relationship between , the refractive structure constant , the wavelength 


of light X, and the distance to target I can be quantified by 


Tq » 0.33 




3/ T \/5 


( 4 ) 


For typical wavelength and target ranges for HEF applications, can range from ~one 
centimeter in strong turbulence to lO’s of centimeters in weak turbulence. The Cl value 
(in units of m'^^^) characterizes the strength of the turbulence. The Cl value was first 
theoretically parameterized by Audrey Kolmogorov and has been measured at numerous 
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locations worldwide [4], The signifleanee of the value stems from Kolmogorov’s 
original theory for turbulent flow, whieh has formed the basis for all eontemporary 
turbulenee theories [4], 

The basie idea behind Kolmogorov’s theory is this: Turbulent fluid flow is 
governed by the Navier-Stokes equations, but signifieant ehallenges he in solving these 
equations within a system as expansive as the atmosphere. Kolmogorov’s theory 
therefore utilizes a statistieal approaeh in deseribing the flow of kinetie energy from 
large-seale eddies, whieh ean be tens to hundreds of meters in size, into smaller eddies, 
whieh are typieally millimeters to eentimeters in size [10]. Eddies within this range of 
sizes, whieh is defined as the inertial subrange, are assumed to be relatively 
homogeneous and isotropie within smaller regions of spaee [4], and thus properties sueh 
as the refraetive index vary over spaee and time in identieally-distributed inerements. 
This assumption allowed for Kolmogorov to utilize a strueture funetion in deseribing 
turbulenee, whieh is represented by the value. 

In atmospherie modeling, an inerease in the value yields an inerease in the 

amount of atmospherie turbulenee. Furthermore, the parameter tends to be a strong 
funetion of altitude, and generally deereases at higher altitudes. Weak atmospherie 
turbulenee would exhibit a value of «10"'’m“^'% while strong turbulenee would 
oeeurif C" «10-'^m“"'\ 

C. THERMAL BLOOMING 

Thermal blooming deseribes a non-linear phenomenon in whieh propagating laser 
light deposits a small portion of its energy along the beam path into the loeal atmosphere. 
As the laser energy is absorbed by the moleeules and aerosols in the air, the ambient air 
temperature inereases, ereating loealized gradients in air density, subsequently altering 
the refraetive index of the medium [11] (see Equation 3). A suffieient ehange in the 
refraetive index effeetively ereates a diverging, eoneave lens within the volume of the 
laser beam, whieh exhibits a greater lens thiekness and thus optieal density at the edges 
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than at the center [5]. Ultimately, this lens causes the beam to “bloom,” or diverge from 
its propagation path, effectively reducing the irradiance on target. 


1. Paraxial Wave Equation 

The effects of thermal blooming on HEL beam propagation can be estimated by 
first assuming that the electric field of the propagating electromagnetic energy takes the 
form of 

E{v) = y/{v)e^^\ (5) 

where is the slowly-varying complex envelope and the harmonic factor has 

been suppressed. Here, k-2nlX is the wave number, X is the laser wavelength, z is 
the primary direction of propagation, and the real part of E specifies the electric field. 
The irradiance of the energy in SI units thus takes the form of 

/(r) = ^£-E = ^rV, (6) 

where y/* represents the complex conjugate of wave function y/ . Equation (5) represents 
a solution to the paraxial wave equation, which follows as [11] 

^ = -^y/ + ^V\y/ + {ik5n)y/, (7) 

where -ey/! 2 represents atmospheric extinction (see section III A.), N]y// 2k 
represents transverse diffraction, / dx^ + / dy^ is the transverse Eaplace 

operator, and the {ik5n)y/ component describes the effects of varying the refractive 
index 5n due to thermal blooming (and in principle, turbulence). 


As stated previously, a change in air temperature has a direct effect on the 
refractive index of the absorbing medium, described by 

■1 


dT T„ 


at. 


( 8 ) 


where is the averaged refractive index of air, is the ambient temperature, and AT 
represents the increase in temperature due to absorption. 
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For most HEL weapon systems, it is common to assume that the rate of heating is 
too small to create pressure differences (assumed to be nearly instantaneous) [11], 
resulting in constant-pressure, or isobaric, heating. The temperature rise due to laser 
heating under these conditions can be described using the energy balance equation, given 
by [11] 


dM 


a 


■/-(vV^)Ar + 


K 


pC 


■v^ 


1 

p J 


M . 


( 9 ) 


dt pC^ 

In Equation (9), the al / pC^ term represents laser heating, where a is the 
atmospheric absorption coefficient, p is the air density, is the constant-pressure 
specific heat, and / is the laser irradiance. The (v»V^) AT term represents the effects of 
energy exchange by convection, where v is the wind velocity, is the transverse 
gradient operator. The / pC^hJ term represents conduction of heat through the 

air, where K is the thermal conductivity. 


In most HEE propagation scenarios, convection dominates conduction except for 
conditions involving absolutely calm winds. If one assumes that the laser beam’s pulse 
repetition frequency (PRE) remains low enough to allow the locally heated air to 
be removed from the beam path by either wind or beam motion, then the conduction 
term can essentially be ignored in what is defined as convection-dominated thermal 
blooming [11]. 


2. Steady-State Phase Distortion without Wind 


In the case of a stationary beam and no wind (|v = O) , the thermal blooming due 


to absorption approaches the steady-state conduction-dominated scenario. The magnitude 
of conduction-dominated thermal blooming on the laser beam can be evaluated using the 
dimensionless scaling parameter given by [12] 


D^=- 


rij-Paz 


^ IzrKna^ ’ 


az <s:l. 


( 10 ) 


17 



where from Equation (8) we obtain that nj.=dn/dT =-[n^I, P is the laser 

power, z is the distanee that the beam has propagated through the absorbing medium, 
and a is the beam radius. The on-axis deerease in beam intensity due to eonduetion- 
dominated thermal blooming ean be approximated by [11]. 

For the eontribution of eonduetion to thermal blooming to be ignored, the wind 
speed |v| must satisfy the condition where ^ = ^1 pC^ represents the 

thermal diffusivity of the medium (z ~ 0.2cm^ / s for air at 300 K) [11]. Because the 
ratio of the thermal diffusivity to the beam radius is virtually negligible in comparison to 
wind speed for all scenarios but those involving absolutely no wind, this condition is 
nearly always satisfied for laser propagation. 

3. Steady-State Phase Distortion with Wind 

Phase distortion due to thermal blooming in the presence of wind has a very 
specific effect on the laser beam irradiance on target, illustrated by Figure 5. 


incident refractive distorted 

INTENSITY TEMPERATURE INDEX INTENSITY 



Figure 5. The effects of steady-state thermal blooming with a transverse wind, 

from [11]. 
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As the transverse wind component moves across the beam propagation path, it 
causes the local air temperature to increase in the downwind direction. As temperature 
change is inversely proportional to both the change in air density and the refractive index, 
both quantities decrease in the same direction as the wind. Because light bends in the 
direction of an increasing refractive index, the laser beam shifts its peak intensity in the 
upwind direction. 


This effect can be roughly quantified by assuming steady-state thermal blooming, 
in which dlsT I dt — Q in Equation (9), and that the wind velocity is entirely along the 

transverse x-direction, so \ -vi. The phase shift (in rad) from thermal blooming can then 
be written as [11] 


tpg[x,y,z) = ik{Sndz = -——- lldx', (11) 

0 PCpVi 

where the beam is assumed to be collimated along the z-direction. If we further assume 
an irradiance distribution for a Gaussian beam. Equation (11) can be solved analytically 
to obtain [11] 




i^y / of l-l-e//(x/a)]. 
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where the maximum phase shift across the beam can be expressed as 




1 


24^ 




and 
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(15) 


19 



The Njy parameter represents the Bradley-Hermann [13] distortion number and is 
utilized as a quantitative measure for the intensity of thermal blooming [11], Blooming 
becomes a significant complication to the propagation of laser energy when > 25 [4], 

The effects of phase distortion on the laser beam can be specifically inferred from 
the series expansion in Equation (13). The leading (x/ a) term produces the deflection of 

the laser beam along the wind direction in the x-direction and the leading (y / term 
produces the beam divergence in the y-direction [11], 


4. Laser Beam Irradiance Pattern Due to Thermal Blooming 


The thermal blooming irradiance pattern can be estimated for the phase 
accumulation by solving the paraxial wave equation in the ray-optics limit using 
perturbation theory. To the first-order perturbation, the solution follows as [11] 


I{x,y,z) 
I e 


= exp 




(16) 


where 1^ is the unperturbed, Gaussian, collimated beam profile and is the collimated 
beam irradiance distortion parameter with the inclusion of wind effects, taking the form 

of[ll] 


, _{nQ-\)aPz^ 
nn^pCTva^ 


(17) 


In general, an increase in the beam distortion parameter represents a greater effect 
of thermal blooming on the laser beam and results in an irradiance pattern taking the form 
of a crescent shape, which is illustrated in Figure 6. 
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Figure 6. The general transverse laser beam profile due to the effects of 
steady-state thermal blooming. The presence of wind moving from left 
to right in the figure results in the typical crescent shape shown. The 
color map indicates the irradiance of the laser beam, with the dark red 
color representing the highest irradiance and the dark blue the lowest. 
The scale of beam deflection depends upon the distance to the target, 
beam power, atmospheric extinction, and other parameters as described 

in the text. 


21 




THIS PAGE INTENTIONALLY LEET BLANK 


22 



IV. ATMOSPHERIC CODE OVERVIEW 


Effective employment of HEL weapon systems requires knowledge of the 
atmospheric optical properties and how these properties affect laser propagation. The 
atmospheric parameters that determine these optical properties are extremely complicated 
and depend on pressure, temperature, humidity, aerosol size composition, etc. 
Eurthermore, all of these parameters are functions of elevation, time of day, and season. 

Due to these complications, it is difficult to fully map out this multi-dimensional 
parameter space for varying altitudes with experimental data. Thus, a means of 
generating optical properties based upon known physical principles for a large range of 
conditions must be employed, which can then be used to evaluate laser beam 
propagation. 

A. MOLECULAR ABSORPTION: 

Among the primary atmospheric optical properties that govern laser beam 
propagation are the four components of the extinction coefficient (See Chapter II, Section 
A). A combination of measurements and quantum theory is utilized to determine the first 
component, the molecular absorption coefficient , which represents the absorption of 
laser energy by molecules in the atmosphere. 

Within the electromagnetic spectrum, the position and strength of absorption 
features at each wavelength arise from quantized molecular energy levels. The 
atmospheric codes which are employed in atmospheric modeling query a database of 
absorption features (such as the HITRAN database), modify the line shapes based upon 
temperature and pressure, and then add them together to get the absorption coefficient at 
a specific wavelength. This method is based on assumptions regarding the natural line 
shape and how temperature and pressure affect the shape. 
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1. Assumptions 

The general form of an absorption ooeffleient line shape is represented by 

( 18 ) 

where A represents the number of atoms per unit volume, S is the line strength, 
/ (v - ) is the normalized line shape at frequeney v , and is the eentral frequeney of 

the line, where /(v-v^) dv -1. One of the most fundamental assumptions regarding 

moleeular absorption modeling is that the natural line shape for eaeh absorption feature is 
based upon modeling a moleeular eleetron as a damped, driven harmonie oseillator, with 
the equation of motion given by [9]; 

F eE ■■ 

— = — = r+ A 7 K'+^r . (19) 

mm 

e e 

In Equation (19), F represents the foree veetor, is the mass of the eleetron, e 

is the eharge of the eleetron, E is the total eleotrie field (ineluding the eontribution from 
the light), f is the eleetron aoeeleration, T^ is the effeetive natural damping eonstant, f 
is the eleetron veloeity, ^ is related to the strength of the restoring foree (whieh is related 
to the effeetive spring eonstant), and f is the eleetron position. While the potential 
energy of moleeular bonds, also known as the Morse potential, is not the same as a 
harmonie oseillator, at lower energy levels it ean be reasonably approximated by the 
harmonie potential, as shown in Figure 7, thus providing motivation for using Equation 
(19). This approximation is nominally satisfied at typieal atmospherie eonditions, as most 
of moleeular eleetrons remain in their ground energy state at standard temperature and 
pressure (STP). 
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Figure 7. Harmonic oscillator approximation. At low energy levels and small 
intemuclear separation, the Morse potential can be adequately 
approximated by the harmonic potential, after [15]. 


Equation (19) can be used to determine the complex refractive index [9]. The 
imaginary part of the refractive index is responsible for absorption and determines the 
functional shape of / (v - ) as 


f{y-yo) = 


1 


( 20 ) 


where T ^ is the same effective damping constant in Equation (19). Equation (20) is in 
the form of a Lorentzian, with E^ also specifying the full-width half-maximum of the 
response, as visualized in Figure 8. 
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Figure 8. The dependence of spectral line width / (v - ) and height on 

resonant frequency and damping, respectively, after [16]. 


2. Absorption Line Broadening 

The absorption spectrum of a single line is broken down into the shape and 
strength of the line as described by Equation (18). However, various mechanisms affect 
these spectral line characteristics, which are discussed below. 


a. Natural Line Width 


The natural line width of an absorption feature is characterized by the damping of 
the harmonic oscillator vibrations due to energy emission of the molecule [9]. The 
Heisenberg uncertainty principle can be used to estimate the width of the absorption line 
and thus provides a means of estimating T,,. It states 


AEAt> — , 

l7t 


Av = 


AE 

~h 


1 

InAt ’ 
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( 22 ) 


where AE is the change in energy of the molecule for a particular transition. At is the 
lifetime of the excited state, h is the Planck constant, and Av is the width of the line 
shape. 
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In the thermal infrared regime, the lifetimes of isolated moleeules are 
At « 0.1 -10 s, in whieh ease the width of the lineshape A v ~ 0.3 Hz, whieh means the 
damping eonstant T^ is similarly small. Beeause the natural line width determined from 

the uncertainty principle and modeled by quantum mechanics is very small, the effects of 
the two broadening mechanisms discussed below dominate the spectral line shapes and 
can be neglected. 


b. Pressure Broadening 


As stated earlier, two broadening mechanisms dominate the characteristic shape 
of spectral lines. Pressure broadening is a mechanism by which molecular collisions in 
the atmosphere transfer energy between molecules, thereby reducing the lifetime of 
higher energy states and increasing the effective damping of the oscillator. It is the 
dominant broadening mechanism in the lower atmosphere below ~ 20 km [14]. 
Increased pressure causes more collisions that broaden the Lorentzian absorption lines, 
taking the form of 


/p(^-^o) 
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where is the effective damping constant due to collisions and is many orders of 


magnitude larger than F^ [14]. The damping constant F^ is determined from the mean 


time between inelastic molecular collisions , so that Fp = 1 / 27i:t^. It also follows the 
trend 
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(24) 


where F^ represents the line width at a reference temperature and pressure , and 
n «1/2 is determined empirically [9]. 


Figure 9, below illustrates the Lorentz line shape for three pressure values. For 
typical vibration-rotation bands, a line width of 1.5x10^ Hz at a pressure of one bar is 
typical [17]. 
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Figure 9. Lorentz line shape eaused by pressure broadening, with variations in 
damping constant and atmospheric pressure, after [14], 


Doppler Broadening 


c. 


Doppler shifts also serve to broaden spectral line shapes. As laser energy 
propagates through the atmosphere, the light encounters molecules with a random 
distribution of velocities. A Doppler frequency shift due to molecular motion results in a 
further broadening of the emission/absorption spectrum, but in a Gaussian profile 
described by [18] 



(25) 



(26) 


where is the Doppler line half-width, is the Boltzmann constant, m is the mass of 
the molecule and c is the speed of light. 

Since a change in temperature implies a change in the velocity of molecular 
motion, an increase in temperature generally correlates to an increase in Doppler 
broadening. 
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d. Voigt Profile 

Atmospheric codes use Voigt profiles to combine the effects of pressure 
(Lorentzian) and Doppler (Gaussian) broadening mechanisms. The Voigt line shape is 
quantitatively described by [14] 




r +00 — \) 

r e ^ 


-dy, 


+r 

where T = and x = (v - Vg) / . The Voigt profile is illustrated in Figure 10. 


(27) 



Figure 10. Collision, Doppler and Voigt spectral line profiles, after [14]. 


In the shoulders of the line shape, the Voigt profile exhibits a Lorentzian shape, while a 
resemblance to the Gaussian shape is seen near the center of the line shape. 

3. Continuum Absorption 

Although the Voigt profile is a well-known representation of spectral line shapes, 
it is not completely devoid of error. Continuum absorption, a phenomenon discovered in 
1918, represents the discrepancy between such data obtained from a mathematical 
summation of Voigt lines and actual measurement data of the molecular absorption 
spectrum. 
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The physical principles behind continuum absorption remain debated; however, a 
couple of theories attempt to reconcile the difference between experimental and modeling 
data. One possible explanation lies in the hypothesis that actual absorption line shapes are 
not perfectly modeled by the Voigt profile, and a summation over thousands of such lines 
gradually yields a rather significant margin of error. Another possible explanation lies in 
the tendency of water to form dimers in the atmosphere. A dimer consists of two H 2 O 

molecules loosely-bound with a hydrogen bond, exhibiting measured dissociation energy 
of ~ 13.2±0.12kJ/mol [19]. Incident laser light upon such dimers could cause the beam 
to deposit extra energy into the atmosphere, resulting in an absorption spectrum which is 
not properly modeled by the Voigt profile. 

Whatever the mechanism, continuum absorption varies slowly with wavelength, 
and thus atmospheric codes often add the contribution from continuum absorption using 
the semi-empirical Clough-Kneizys-Davies (CKD) continuum model [20]. 

B. MIE SCATTERING THEORY: , /?,„, AND /?„ 

While a combination of measurements and quantum theory is utilized to 
determine the first component of the extinction coefficient, a separate solution called Mie 
scattering is used to determine the latter three components. 

1. Theory Introduction 

Mie scattering theory describes the interaction between electromagnetic radiation 
and small particles suspended in the atmosphere, such as dust and aerosols [9]. The 
theory is utilized to determine the molecular scattering coefficient as well as the 

aerosol absorption coefficient and aerosol scattering coefficient . Mie theory 

contains three key assumptions, the first of which states that all atmospheric particles are 
spherical in shape. The second assumption states that the scattering particles are 
homogeneous in material and thus can be represented by a single complex refractive 
index h. The last assumption holds that the light incident upon a particle, in this case 
from the laser, is in the form of a plane wave. 
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While these assumptions simplify the interaetion of light and atmospherie 
partieles, they do not always represent reality. For example, many partieles sueh as eubie 
salt erystals are not spherically shaped. However, Mie scattering theory can still be 
applied in that case since the irregular particles are suspended in a turbulent media and 
are not preferentially oriented in one direction [21]. More sophisticated codes extend Mie 
scattering to include spheres with layers of thin films for various purposes, such as 
modeling water layers on aerosol particles [22]. 

2. Amplitude Scattering Matrix 

In calculating the atmospheric extinction coefficients, the initial goal of Mie 
scattering theory is to determine the scattered field amplitude at an angle 0 and a 
specified distance r from a scattering object, in this case a suspended particle or aerosol 
[9]. While Mie theory utilizes Maxwell’s equations to determine the solutions to the 
fields, both inside and outside of the scattering object or scatterer, HEL propagation only 
concerns “far-field” scattering, where r » / /I and a is the particle radius. In this case, 

the fields take the form of [9] 


^e,(r,0) 

exp(-/A:r) 

's^{e) 0 
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where E''] and Ef are the parallel and perpendicular components of the incident electric 

field and and E^ are the respective components of the scattered electric field. The 
wavenumber k-lnlX is representative of light in air since the refractive index of air 


h — \. The amplitude scattering matrix 


's^{e) 0 

0 


determines the amount and 


polarization of scattered light at a given angle 6 measured with respect to the incident 
light. 


The scattering amplitudes 5j and take the form of [23] 
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(30) 


5,(e)=Z 

l=\ 

where n^ and depend only on the scattering angle 0 and are defined as [23] 
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In Equations (31) and (32), is the associated Legendre polynomial. The Mie 


coefficients and are dependent only on the optical wavelength A, the scatterer 


radius a , and the scatterer index of refraction n , and are defined as [23] 
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(34) 


where ^and denote the Hankel functions [19]. 


3. Total Cross Sections 

The total scattering cross section integrated over all angles takes the form of [9] 

n 2 cc 

<^.=5^Z(2< + 1)(KI +|(’,l )■ (35) 

The extinction (scattering + absorption) cross section is given by 

2 CO 

<7,=—X(2( + l)SH(a,+(.,), (36) 

LTl 

where 9^(...) denotes the real part. Finally, the total absorption cross section is [9] 

(37) 

which is simply the difference between the extinction and scattering cross sections. The 
absorption and extinction coefficients are calculated from their respective cross-sections 
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by integrating over the contributions from different particle sizes. For example, the 
aerosol scattering coefficient [5^ is calculated by integrating [9] 

Pa{^) = \lN{r)c7^{r)dr, (38) 

where N dr is the number of particles with a radius between r and r + dr per unit 
volume. 


4. Rayleigh Scattering Limit for Determining 


Molecular scattering assumes that particles of size a are much smaller than the 
wavelength of the laser light In the Rayleigh scattering limit, the scattering 

cross-section depends on the wavelength of electromagnetic energy as qc , so shorter 
wavelengths scatter more readily. A famous example of this principle is the color of the 
daytime sky: because blue light is of a shorter wavelength than red light, it Rayleigh 
scatters more readily and thus the sky appears blue. 


In this case, the Rayleigh scattering approximation can be utilized to determine 
the scattering cross section by [9] 


128;r' a 
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(39) 

(40) 


The molecular scattering coefficient can then be calculated from the cross section 
following the same method as Equation (38). 


Figure 11 visualizes the general scattering cross section-dependent trends 
associated with Mie scattering theory. In the low-frequency regime the circumference of 
an arbitrary spherical particle is smaller than the wavelength of energy passing over it, 
and thus Rayleigh scattering ( oc ) describes the scattering cross section . 
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Figure 11. Mie theory scattering behavior with respect to the cross-section of 
the scatterer. The Rayleigh scattering limit is indicated in the figure. 

5. Aerosol Absorption and Scattering Coefficients and 

Aerosol sizes are typically on the order of an optical wavelength (a «A). 
According to Mie theory, aerosol-sized particles scatter and absorb electromagnetic 
energy in a manner weakly dependent on the wavelength of light passing through them 
[9]. Figure 11 visualizes this trend; In the high frequency (optical) limit characterized by 
the resonance when a~ A, the cross section is roughly independent of the wavelength. 
As previously stated, the aerosol absorption and scattering coefficients are calculated 
following the method in Equation (38). Flowever, unlike the Rayleigh scattering 
approximation for calculating aerosol coefficients have no similar approximating 
solutions. 

C. MODTRAN 

There are several well-established radiative transfer codes that are used to 
generate extinction coefficients needed to model laser propagation. One of the most 
widely used of these is MODTRAN (MODerate resolution atmospheric TRANsmission) 
[25], a band-transmission code, which utilizes a database of pre-calculated coefficients 
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for a wide variety of atmospherie eonditions and aerosol eontent. MODTRAN uses the 
HITRAN database for moleeular absorption speetra [26]. 

MODTRAN is a relatively fast (but potentially less aeeurate) method of obtaining 
optieal properties eompared to other mueh more eomputationally intensive eodes sueh as 
the line-by-line LBLRTM program [27]. The reason line-by-line eodes are more CPU 
intensive is that they first sum over all absorption lines, then add eontinuum absorption 
for a speeifie wavelength, then ealeulate Mie absorption and seattering eoeffioients from 
physieal prineiples to finally determine the extinetion eoeffieient. 

However, while MODTRAN may be a potentially less aeeurate means of 
generating optieal properties, the program is very flexible as it allows for various 
atmospherie parameters to be set by the user. 

1. Preset Atmospheric Models 

Table 2 lists the six atmospherie model presets that MODTRAN ean use to 
generate atmospherie optieal properties. The user ean utilize one of these presets as is, 
modify aspeets of these presets, or speeify user-defined models. 


Atmospheric Model # 

Description 

1 

Tropical Atmosphere (15° Latitude) 

2 

Mid-Latitude Summer (45° North Latitude) 

3 

Mid-Latitude Winter (45° North Latitude) 

4 

Sub-Arctic Summer (60° North Latitude) 

5 

Sub-Arctic Winter (60° North Latitude) 

6 

1976 US Standard Atmosphere 


Table 2. MODTRAN preset atmospherie models, after [25]. 
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The six preset models vary in several eharaeteristies, most importantly in their 
temperature, HjO, and O 3 profdes. Figure 12 displays the temperature profdes as 
funetions of altitude for eaeh of these presets [25]. 



Figure 12. MODTRAN model atmosphere temperature profiles. The varying 
eolors indieate temperature profiles for the different model 
atmospheres, from [25]. 


Naturally, the Sub-Aretie Winter preset exhibits the eoldest surfaee temperature 
while the Tropieal model exhibits the warmest. 

Figure 13 displays the HjO and O 3 profiles as funetions of altitude for the six 
preset atmospheres [25]. 
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Figure 13. MODTRAN model atmosphere H 2 O and O 3 profiles. The gaps in the 
data are the result of the breaks in the vertieal and horizontal axes, as 

indieated, from [25]. 


As expeeted, the Tropieal model atmosphere eontains the highest eoneentration of H 2 O 
at the surfaee and the Sub-Aretie Winter eontains the lowest. The O 3 profile is slightly 
more eomplieated, with all models eontaining a roughly equal eoneentration at ground 
level, and varying greatly at the peak at ~ 20 km, whieh represents the ozone layer. 


2. Meteorological Range Presets 


MODTRAN allows for varying aerosol densities by speeifying the ground 
meteorologieal range presets at 23, 15, 10 or 5 km. Meteorologieal range (MR), in this 
ease, is defined by Equation (49) as [25] 

MR[km] =-(41) 

EXT550[km-'] + 0.01159km-‘ 

where EXT550[km''] represents the surfaee aerosol extinetion at 550 nm and 


0.01159km"‘ is the surfaee Rayleigh seattering eoeffieient at 550 nm. Meteorologieal 
range ean also be equated to the more subjeetive visibility, defined as the maximum 
distanee at whieh an individual ean spot and identify a prominent objeet along the 
horizon at day (or a light at night), by the approximate Equation (42) 


Visibility « 


MR 

1.3T0.3 ■ 


(42) 
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If no meteorological range is specified, a default value is generated from presets 
associated with the boundary layer chosen in the particular aerosol model (see next 
section). 


3. Aerosol Models 

Another option within MODTRAN specifies the aerosol model preset used for the 
surface layer (from zero to two kilometers above sea level) that also includes a default 
surface-meteorological range. While MODTRAN also allows the user to specify aerosol 
models at altitudes above the surface layer, this thesis will only consider engagements 
below two kilometers. The surface layer atmospheric extinction coefficient dependence 
on the relative humidity is determined by the water-vapor content of the selected 
atmosphere model. Table 3 lists the primary aerosol presets available in MODTRAN, 
with model #7 as the user-specified aerosol extinction coefficient option excluded. 


Aerosol 

Model # 

Description 

0 

No aerosol or cloud attenuation 

1 

RURAL model; Default meteorological range set to 23 km 

2 

RURAL model; Default meteorological range set to 5 km 

3 

NAVY MARITIME model; Default meteorological range set based on wind speed & relative 

humidity 

4 

MARITIME model; Default meteorological range set to 23 km 

5 

URBAN model; Default meteorological range set to 5 km 

6 

TROPOSPHERIC model; Default meteorological range set to 50 km 

8 

FOGl model (Advective Fog); Default meteorological range set to 0.2 km 

9 

FOG2 model (Radiative Fog); Default meteorological range set to 0.5 km 

10 

DESERT model; Default meteorological range set by wind speed 


Table 3. MODTRAN main aerosol and cloud options, after [25]. 


38 














The preset most commonly used in this study was the #4 Maritime model, with a 
default ground meteorological range of 23 km. The primary purpose of the Maritime 
model is to study the propagation of infrared energy through the atmosphere in the near 
sea surface layer, an area in which all Naval vessels would operate in, making the model 
useful for this research. 

This model utilizes an aerosol profile composed of a sea-salt component and a 
continental component, the latter of which is primarily comprised of water soluble 
material such as ammonium, calcium sulphate, and other organic compounds. This 
aerosol profile does not contain any larger dust particles due to the assumption that these 
particles would rapidly fall out of the atmosphere as a result of gravitational settling, 
which occurs on the order of hours to days as the air moves from land and across the 
ocean [21]. 

Furthermore, the Maritime model also neglects the effects of fresh sea-spray, 
which occurs in the lower altitudes above the sea surface at 10-20 meters and is strongly 
dependent on surface wind speed. Other Maritime aerosol models are available in 
MODTRAN that include these sea-spray effects and could represent a possible extension 
to this research. 

4. Cloud and Rain Models 

Additional options specify the cloud and rain presets utilized by MODTRAN in 
generating atmospheric extinction coefficients. The rain profiles decrease linearly from 
the ground to the top of the associated cloud model, with the rain cutting off at the cloud 
top [20]. Table 4 lists the major MODTRAN cloud and rain models, with preset # 11 as 
the user-specified cloud extinction and absorption coefficient option excluded. 
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Rain Model # 

Description 

0 

No clouds or rain 

1 

Cumulus cloud layer with base at 0.66 km and top at 3.0 km 

2 

Altostratus cloud layer with base at 2.4 km and top at 3.0 km 

3 

Stratus cloud layer with base at 0.33 km and top at 1.0 km 

4 

Stratus/stratocumulus cloud layer with base at 0.66 km and top at 2.0 km 

5 

Nimbostratus cloud layer with base at 0.16 km and top at 0.66 km 

6 

2.0 mm/hr ground Drizzle (Includes ICLD # 3 with 0.86 mm/hr rain at 1.0 km) 

7 

5.0 mm/hr Light rain (Includes ICLD # 5 with 2.6 mm/hr rain at 0.66 km) 

8 

12.5 mm/hr ground Moderate rain (Includes ICLD # 5 with 6.0 mm/hr rain at 0.66 km) 

9 

25.0 mm/hr ground Heavy rain (Includes ICLD # 1 with 0.2 mm/hr rain at 3.0 km) 

10 

75.0 mm/hr ground Extreme rain (Includes ICLD # 1 with 1.0 mm/hr rain at 3.0 km) 

18 

Standard Cirrus model with 64 pm mode radius for ice particles 

19 

Sub-visual Cirrus model with 4 pm mode radius for ice particles 


Table 4. MODTRAN cloud and rain models, after [25]. 


D. LEEDR 

LEEDR (Laser Environmental Effects Definition and Reference) is a MATLAB- 
based atmospheric modeling program created by the Air Force Institute of Technology 
(AFIT), and serves a similar function as MODTRAN in characterizing a location’s 
atmosphere based on various user inputs such as location, weather, etc. [28]. Because 
using MODTRAN to generate optical parameters for model atmospheres containing 
variable amounts of precipitation proved unsuccessful, LEEDR was thus implemented to 
obtain this data set. All LEEDR-generated rainy model atmospheres were calculated 
using atmosphere types and aerosol models corresponding to those available in 
MODTRAN. While there were slight variations in the final propagation results obtained 
from the model atmospheres generated by MODTRAN and LEEDR, the data strongly 
agreed at the 1.064 pm wavelength (which corresponds to the nominal laser wavelength 
for the SSL). 
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V. SIMULATION METHODS 


In order to accurately and efficiently model HEL propagation in various 
atmospheric conditions, NFS has developed ANCHOR (Atmospheric NFS Code for HEL 
Optical pRopagation). ANCHOR is a succinct Matlab script that utilizes pre-calculated 
data sets containing extinction coefficient values across a specified wavelength spectrum. 
The components of the extinction coefficient - namely the molecular and aerosol 
absorption and scattering coefficients (see Chapter III, Section A) - and specified laser 
parameters are subsequently integrated into the scaling equations governing laser 
propagation. The on-target laser irradiance and power-in-the-bucket values are obtained 
across a specified wavelength spectrum. 


Eour specific effects govern the behavior of the beam as it propagates to the 
target: Diffraction, turbulence, platform jitter, and thermal blooming. The linear effects of 
diffraction, turbulence and jitter are assumed to produce a time-averaged Gaussian 
irradiance profile, and can be evaluated independently of each other. These effects are all 
implemented within the “master Equation (43) [29],” which describes the time-averaged 
irradiance, or intensity, (/) of the laser beam on a target (in W/m^). 


n{w. 


2 ^ ^TB • 


(43) 


In Equation (43), represents the total output power from the laser beam 


director, is the time-averaged radius (1/e in the intensity) of the laser spot on 

target. The factor comes from Beer’s law (see Chapter III, Section A), where e is the 

total extinction coefficient due to atmospheric absorption and scattering (in m"') and f is 
the distance (in m) from the beam director to the target. The factor represents the 

thermal blooming Strehl ratio, which accounts for the non-linear reduction in the on-axis 
irradiance due to thermal blooming. 


The total laser beam spot size on target is determined by [29] 


( 44 ) 
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in which represents the root-mean-square (rms) spot size due to only diffraction in a 

vacuum, the rms spot size due to only turbulenee, and w- the rms spot size due to 

only platform jitter. Sinee these effeets are independently assumed to result in time 
averaged Gaussian irradianee patterns, their eombined rms eontributions add in 
quadrature. 

A. LINEAR EFFECTS 

Before the linear effeets of diffraetion, turbulenee and jitter governing the 
behavior of a propagating laser beam can be implemented into Equation (43), they have 
to be first eonsidered and quantified independently. 


1 . 


Vacuum Diffraction: 


The effeets of laser beam diffraetion in a vaeuum can be quantified by estimating 
the beam’s spot size due to diffraetion, whieh follows as 

lU 


w. ~ M- 
“ kD 


(45) 


where D represents the beam diameter at the souree and M represents the mode faetor 
measuring beam quality. If the laser has a uniform transverse beam shape, then D is also 
the diameter of the beam direetor; for a Gaussian shape, it represents the diameter of the 
beam at the He point in the irradianee. The mode faetor is a unit-less parameter whieh 
represents the mode eontent of a multimode laser beam obtained from a stable resonator. 
The square of M, the parameter, compares the produet of beam divergenee 0 (in 
rad) and spot size >v of a real laser beam to that of a fundamental Gaussian referenee 
beam, so that [30] 

0w 


=■ 


0 , 

ref ref 


(46) 


where 0^^^^ represents the referenee beam divergence and the referenee beam spot 

size. The parameter is a widely used quality metric for laser beams, and as sueh is 
often defined as the “beam quality” faetor. However, a more aeeurate definition of is 
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the “beam propagation” faetor, since the parameter ultimately provides a measure of how 
the beam will propagate through free space [31]. For an ideal Gaussian beam, =1, 

while for a non-ideal beam, > 1. An value less than unity is not possible. 


The spot size due to diffraction can be further described by taking into account the 
laser beam’s focal range F , which yields 

v 2 


2 

=-^ 2—2 


k^w 


+ w 11 - 


(47) 


where the wavenumber k -Itt / A and represents the laser beam radius at the source. 


For this research, the focal point of the laser beam was placed on the target, so F = £ and 
the second term in (47) is eliminated. 


2. Atmospheric Turbulence: w, 


Turbulence is described comprehensively in Chapter III, Section B. The effect of 
turbulence on the laser spot size on target can be estimated by 

Ini 






(48) 


where I once again represents the distance to target and is the Fried parameter, which 

defines a length scale within which the laser beam maintains transverse coherence during 
its propagation path. The Fried parameter is estimated by 


r„«3.03 


V' 


/ly 


- 0.6 


C^l 


(49) 


where characterizes the strength of the turbulence, as discussed in Chapter III, 
Section B. 


3. Platform Jitter: Wj 

Platform jitter is defined as the randomized movement of the laser spot position 
on the target due to the motion of the laser platform as well as tracking error by the active 
or passive tracking system. The time-averaged effect of jitter can be quantified by 
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considering Equation (50), where 8^^^ defines the total rms angular displaeement (in rad) 
due to all eontributions from jitter and t is the distanee (in m) to target; 


w. 

} rms 


8 L 


(50) 


B. THERMAL BLOOMING 


As previously discussed, thermal blooming deseribes the non-linear phenomenon 
in whieh ambient air surrounding the propagating laser beam is heated by the beam itself, 
resulting in a ehange in the loeal air density. This ehange in air density results in a ehange 
in the refraetive index of the loeal atmosphere, ereating a “lens” in the air that eauses the 
beam to diverge. The direet effeet of thermal blooming is the reduetion of the laser beam 
irradianee on target. Thermal blooming is diseussed extensively in Chapter III, Seetion C. 


1 . 


Distortion Number: N, 


The effeets of thermal blooming can be charaeterized by the thermal distortion 
number (Chapter III, Seetion C.3), whieh is defined as [4]; 





r«l 

( \ 

ir| 

(^) 


(^) 

Pf, \ 

1 v| 



(^) 



dz. 


(51) 


where the wavenumber k-lnlX, ~ 1 kg/m^ represents the ambient air density, 

Cp ~10^ J/kg-K is the speeifie heat of air at eonstant pressure of air, cc(z) is the 

absorption coeffieient (in m"'), r(z) = represents the total transmission of the laser 

energy, rij. (^z) = dn/ dT « -10^® K ' is the change in the refraetive index with respeet to 

temperature, v(z) is the wind speed in the transverse direetion of the beam, and D(z) is 

the beam diameter. In all of the above expressions, z is the distanee from the beam 
direetor and is measured along the propagation path. 


a. Constant Atmospheric Properties 

Given a horizontal beam path from the souree to the target, whieh is valid only if 
both the souree and target are at equal altitudes above the sea level, the atmospheric 
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absorption coefficient, refractive index, and wind speed are all assumed constant along 
the propagation path, and thus the distortion number can be approximated as 


b. 



Beam Diameter: Di^z) 


(52) 


The beam diameter Di^z) along the propagation path has to be estimated in order 

to evaluate Equation (51). At the beam director, the laser beam is assumed to have a 
diameter of . At the target (assumed to be at the beam waist), the beam has a time- 

averaged radius defined by Equation (44). Assuming that the beam propagates between 
these two points and the fundamental mode is roughly Gaussian, the beam diameter along 
the path can be estimated as 

D(z)^2w ll + zY, . (53) 

V / 


In (53), represents the laser beam spot radius at the target, z is the distance from the 
target and is the “effective” Rayleigh range, which is defined such that 

E)(z = -f) = f)o [5]^ and visualized in Eigure 14. Eollowing this definition, the 
“effective” Rayleigh range can be described by 
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Figure 14. 


Visualization of laser beam profile when Z)(z 


^0 ^ . After [6]. 


-£) = Dq and 


However, if , then the waist of the beam is located at the beam director 

(now at z = 0) and the target is located at z = i, which is visualized by Figure (15). The 
beam diameter is then calculated as 


where 


D{z) = D, 



^eff 



if D,<2w,^,. 


(55) 


(56) 


X 



Figure 15. Visualization of laser beam profile when target is located at z = f 

and Dq < , after [6]. 
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2 . 


Strehl Ratios: S 


The thermal blooming Strehl ratio Sjg represents a ratio of the peak irradianee of 
a laser beam to the peak irradianee of a referenee beam where no thermal blooming has 
oeeurred. This ratio is utilized by the ANCHOR seript to ineorporate thermal blooming 
effeets into the time-averaged irradianee and power-in-the-bueket ealeulations [29]. The 
Strehl ratio for a Gaussian laser beam shape can be estimated by [4] 

Sg= -^-r, (57) 

1 + 0.0625A^ 

while the Strehl ratio for a uniform, or “top-hat,” laser beam shape can be calculated from 


1 + O.OlA)," 

where once more represents the thermal distortion number defined by Equation (51). 
For the propagation simulations in this research, the beam shape described by (58) was 
utilized. 


C. IRRADIANCE LINEWIDTH CONVOLUTION 

While Equation (43) adequately calculates time-averaged irradianee profiles for 
various laser beam configurations over a broad wavelength spectrum, the calculation by 
itself does not take into account the finite linewidth of a laser (—0.1-1.0% for HEEs). 
This characteristic is visualized in Figure 16, in which the horizontal axis represents a 
small portion of the electromagnetic spectrum and the vertical axis defines the strength of 
molecular absorption for each wavelength. The red line indicates a typical 1% HEE 
linewidth. As evidenced by the figure, many detailed absorption features are contained 
within the linewidth, and the overall on-target irradianee profile of the beam is greatly 
dependent on how these features affect the laser. 
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Wavelength (pm) 

Figure 16. Absorption features contained within a finite 1% laser linewidth. The 
red line indicates the spectral range of the HEL, after [26]. 

In order to more realistically model laser irradiance profiles, the variations in the 
irradiance due to the absorption features within a finite laser linewidth must be averaged. 
Considering a distribution of calculated irradiance profiles across a specified frequency 
spectrum, a typical linewidth can be described by the change in laser power as a function 
of the change in frequency dP / df, in which 



Combining Equations (59) and (43), the general formula for the time-averaged irradiance 
for an arbitrary linewidth becomes 

= (60) 

tot 

To simplify the above equation, the assumption is made that the laser frequency spectrum 
has a uniform (“top-hat”) distribution centered at frequency and with a linewidth A/ , 

as shown in Figure 17. This simple top-hat approximation can be improved upon in the 
future, but for now adequately captures important spectral features. 
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Figure 17. Uniform distribution of a laser line eentered at frequeney with 

linewidth A/ . 


Following Equation (43), the time-averaged and linewidth-averaged irradianee 


then beeomes 



(61) 


In implementing this ealeulation, whieh amounts to a “moving average” of irradianee 
values aeross a full frequeney speetrum, Equation (61) is diseretized in order to yield 



(62) 


where Sf is the frequency difference between adjacent sampling points and A/// 
represents the aforementioned linewidth. Typically the linewidth is ~ 0.1-1.0% for 
HELs. 

D. ADDITIONAL HEL PERFORMANCE METRICS 

While irradianee on target is an important parameter to evaluate, it is not the only 
way to estimate laser performance and damage. 
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1 . 


Power-in-the-bucket: Pg 


Power-in-the-bucket is defined as the amount of power from the HEL that 

falls within a pre-defined circle (the “bucket”) with radius on the target. Thermal 
blooming causes the laser beam to smear out in a complex manner, so an effective beam 
size w,ot^eff is estimated as 


(63) 


w. 


tot,eff 


W., 


where w,", is determined from linear effects. The irradiance can then be written as 



-ft 


nw. 


tot,ejf 


(64) 


If , then the power-in-the-bucket is approximated by 

= ( 65 ) 

which is to say that all the power reaching the target plane falls within the bucket. If 
h < ^tot,ejf ’ power-in-the-bucket can be simply approximated by 

«(/)(;rr/). (66) 

2. Dwell Time: 

A more physical parameter that can be used to evaluate laser beam performance is 
the dwell time for melting a through a certain amount of target material. Given an 

arbitrary material of a specified thickness, one could estimate the length of time 
necessary for a laser beam to remain on target in order to melt that given volume of the 
material. 


a. Conductive and Radiative Power Transfer/Loss 

As a laser beam dwells on a target, a certain amount of the laser’s energy is lost 
due to energy conduction away from the target into the surrounding material. This 
conductive power transfer loss (in W) can be quantified by [6] 
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^^ J.rdK(T,-T„) ^ ( 67 ) 

in which r represents the radius of the on-target laser spot (in m), d is the thickness of 
the material, k is the thermal eonduetivity of the material, represents the melting 
temperature of the target, Tq represents the ambient and initial temperature, and L 
represents the distanee through the material over whieh the temperature transitions from 
to Tq . This distanee depends on the material properties and the dwell time; for typieal 

metals and dwell times of ~ 10 s, L~ 2 em [6]. 


The second power loss meehanism is radiative power loss , whieh oeeurs when 

the heated target area radiates as a black body. The radiative power loss can be estimated 
by [6] 

P,=crsA{T:-T:), (68) 

where <7 is the Stefan-Boltzmann eonstant (-5.670x10^* W m"^ K "*), s is the 
emissivity of the material (dimensionless), and the radiating surfaee area A = . 


The total power loss is thus defined as [6] 


Pl- Pc+P r- 


(69) 


In order to caleulate the total power loss , the spot size r must first be estimated. From 
Equation (63), the effeetive laser spot size is estimated as 


r = w. 






\J^TB 


(70) 


This radius depends on the wavelength of the laser, atmospherie eonditions and the 
platform/target geometry. 


b. Time Required to Melt the Target 

Utilizing the previous eomputations, the usable melting power delivered to the 
target ean be quantified by 

P„,=P„-Pl. ( 71 ) 
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in which represents the laser power delivered to the target. Finally, if > 0, the 


dwell time can be quantified by [6] 




^D=- 


melt 


(72) 


in whieh represents the energy required to melt a given amount of the target. 


The dwell time depends on the total range to the target since the melting power 
delivered to the target decreases as the range to the target increases. One can further infer 
that if the melting power delivered to the target does not overeome the power loss due to 
conduction, the dwell time effeetively approaches infinity, in which case the laser beam 
will never significantly heat the target material to its melting point. 


3. Discussion 

To evaluate the effectiveness of laser beam damage to a target, neither the beam 
irradianee nor the power-in-the-bueket ean be independently eonsidered. In fact, both are 
helpful to fully charaeterize the damage. The time-averaged irradianee provides valuable 
information beeause it quantifies the maximum intensity of a laser beam on a speeified 
target and is irrespective of the size of the target. However, the irradianee does not 
indicate by itself how much power is being delivered to the target area. 

Conversely, the power-in-the-bueket data provides a measure of how mueh power 
emitted by a laser beam is eonfined to a partieular spot on a target. Beeause a substantial 
spot size on target is required in order to inflict significant damage, Pg provides valuable 

insight into the lethality of the laser beam. However, this metric depends on the size of 
the bueket selected, under consideration. 

Dwell time, which depends on the material properties of the target and the melting 
volume needed for a kill, can provide guidance for operational procedures in various 
engagement scenarios. Together, all three metrics complement each other and provide a 
well-rounded characterization of laser performance. 
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E. MODEL VALIDATION 


In order to ensure that the ANCHOR code produces realistic laser propagation 
results, irradiance values at specific wavelengths were compared to the irradiance results 
obtained from Wavetrain. Wavetrain is a propagation code produced by MZA Associates 
Corporation based in Albuquerque, NM and Dayton, OH, and differs from the ANCHOR 
script as it actually simulates laser beam propagation by solving the paraxial wave 
equation [32], This method is self-consistent and potentially more accurate than the 
scaling method used in ANCHOR, but also several orders of magnitude slower. 

Wavetrain allows for irradiance calculations to be carried out in model 
atmospheres with thermal blooming effects turned on or off, and similar data can be 
obtained from the ANCHOR script. Irradiance values from Wavetrain were determined 
by averaging its output over many time steps after the pattern reached steady-state, which 
was ~ 0.1 s. 

Table 5 displays ANCHOR and Wavetrain time-averaged irradiance profiles at 
various atmospheric turbulence values taking into account the effects of turbulence and 
thermal blooming for the /I = 1.064 pm wavelength, which corresponds to the standard 
wavelength for a SSL. 
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Value 

Data Type 

Wavetrain (W/m^) 

ANCHOR 

(W/m^) 

0 

Without TB 

3.6x10^ 

3.6x10* 

With TB 

1.7x10** 

1.4x10* 

10-16 

Without TB 

2.5x10* 

1.9x10* 

With TB 

1.3x10* 

8.3x10’ 

lO'® 

Without TB 

6.9x10^ 

2.5x10’ 

With TB 

4.1x10^ 

1.3x10’ 

lO*'' 

Without TB 

4.1x10^ 

1.7x10® 

With TB 

4.9x10® 

1.2x10® 


Table 5. Irradiance propagation data comparison for 1.064 pm laser 

beam. 


Table 6 displays the identical parameters to Table 5, but for the A = 3.82 pm 
wavelength. 


C’ Value 

Data Type 

Wavetrain (W/m^) 

ANCHOR 

(W/m’) 

0 

Without TB 

3.5x10’ 

3.5x10’ 

With TB 

1.9x10’ 

1.9x10’ 

10-16 

Without TB 

3.4x10’ 

3.4x10’ 

With TB 

1.9x10’ 

1.8x10’ 

10 

Without TB 

2.6x10’ 

2.2x10’ 

With TB 

1.5x10’ 

1.2x10’ 

10'^ 

Without TB 

6.0x10® 

3.2x10® 

With TB 

5.8x10® 

2.2x10® 


Table 6. Irradiance propagation data comparison for 3.82 pm laser 

beam. 
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As evidenced by the compared data sets, at small turbulence levels the ANCHOR 
and Wavetrain results exhibit negligible deviation. As turbulence increases to 
« 10"'"'m"^'^, ANCHOR and Wavetrain differ by a factor of ~ two, with the irradiance 

predicted by ANCHOR consistently smaller between the two codes. On closer inspection, 
it appears ANCHOR is overestimating the degradation of intensity due to turbulence 
relative to Wavetrain. Even so, despite the differences that arise between the codes under 
these conditions, both methods predict laser irradiance values well within an order of 
magnitude of each other. It should be noted that Wavetrain takes ~ one minute of 
computation time per run; in that same time, ACHOR completes ~ 10^ runs. 

A trend in the data arises from the observation that as the strength of turbulence 
increases, the strength of thermal blooming decreases. This is due to turbulence in the 
atmosphere causing the laser beam to spread out, lowering the peak irradiance of the 
beam and preventing thermal blooming from taking effect. 
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VI. RESULTS 


The simulated laser propagation results obtained from ANCHOR were analyzed 
and eatalogued in various plots. Due to the effieiency of the ANCHOR seript, each 
generated data plot represents tens to hundreds of thousands of simulation runs; nearly 10 
million runs in total were conducted for this thesis. 

A. STANDARD MODEL PARAMETERIZATION 

To effectively utilize the ANCHOR code in evaluating laser performance, 
standard laser and atmosphere models were chosen in order to establish a baseline 
performance profile for the various studies conducted. Table 7 displays the laser and 
target parameters that define the standard model utilized in the ANCHOR simulations. 


Laser Output Power P 

100 kW 

Laser Wavelength X 

0.5 pm - 5 pm 

Laser Linewidth ^'/ f 

1% 

On-Target Bucket Size 

5 cm 

Beam Quality M' 

3 

Beam Diameter at Source D 

0.5 m 

Platform Jitter 6 

rms 

2.0x10^ rad 

Target Range 

5000 m 

Engagement Geometry 

Horizontal; 10 m Above Water Surface 

Adaptive Optics 

None 


Table 7. Standard model laser and target parameters for ANCHOR. 
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Table 8 displays the atmospherie parameters that define the standard model 
utilized in the ANCHOR simulations. 


Atmosphere Type 

Tropical, 15“ Latitude [MODTRAN Model #1] 

Aerosol Composition 

Maritime [MODTRAN Model #4] 

Temperature, Tq 

300 K (80“ F) 

Absolute Humidity 

19.7 g/m^ (80% Relative Humidity) 

Meteorological Range 

23 km (Clear Day) 

Precipitation 

No Rain 

Turbulence C“ 

n 

(Moderate) 

Wind 

Transverse, 10 m/s 


Table 8. Standard model atmospherie parameters for ANCHOR. 


Finally, Figure 18 displays a simple visualization of the standard model target 
engagement seenario, in whieh both the laser platform and target are positioned 10 m 
above the water surface and are separated by a range of 5000 m. Once again, the diameter 
of the laser beam at the source (and thus the beam director diameter) is 0.5 m. 



Figure 18. Visualization of standard model target engagement scenario. 
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B. STANDARD MODEL STUDIES AT 5 KM TARGET RANGE 


Evaluations of the standard model were conducted by first selecting a single 
parameter (such as turbulence, haze, output power, etc.) and varying it over a realistic 
range while keeping the other parameters fixed to the standard model. The associated 
irradiance and power-in-the-bucket profiles were then determined for each simulation. 
This method allowed various known trends associated with laser propagation to be 
verified and visualized. 

In order to ground the visualized laser performance to realistic operational 
requirements, a “melt threshold” is shown in the irradiance plots. This melt threshold 
represents the minimum on-target laser irradiance required to melt a certain volume of 
the target material, and is a function of the desired area to damage on the target and the 
laser dwell time. In this case, the target was assumed to be an arbitrary aircraft 
(unmanned air vehicle, helicopter, etc.) with a body comprised of 1 mm-thick aluminum. 
In order to attain the melt threshold, the laser must deliver enough energy to the target for 
the aluminum to reach its melting temperature and continue to deliver enough energy to 
melt through the skin, all within a 10-second laser dwell time. 


Table 9 displays the various properties of aluminum necessary for determining 
this melt threshold. 


Density p 

2700 kg/m^ 

Melting Temperature 

933 K 

Heat Capacity 

897 J/kg-K 

Latent Heat of Fusion H 

400 kJ/kg 

Thermal Conductivity k 

237 W/m-K 

Emissivlty a 

0.05 


Table 9. Physical properties of aluminum, from [33]. 
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Estimating the melt threshold requires knowledge of the target material mass m 
and the total energy neeessary to melt this mass , whieh is eomprised of the energy 

required to reaeh the material’s melting temperature Qteaf and the energy required to melt 
the target onee the melting temperature is reaehed . Assuming the desired region to 

melt is eireular with a radius r - 0.05 m and material thiekness d - 0.001m, the mass m 
of the targeted material ean be estimated by 

m = pnr^d = 0.021 kg. (73) 

Next, the energy required to heat the material is [6] 

e„„=mC,(J:.-r.) = 12.0kJ, (74) 

where = 300 K is the ambient temperature. The energy required to melt the aluminum 
onee the melting temperature has been reaehed is 

G...=m/7 = 8.4kJ. (75) 

Finally, the total energy needed to melt is the sum of the two energy values, 

e.,=a..+e„..=20.4kJ. (76) 

With knowledge of the total energy requirement, the on-target laser irradianee 
required to meet the melt threshold ean be estimated by [6] 

(0 ^ 

—+ (Pc+^’,) 

^- - = 0-56 MW/m", (77) 

nr 

where is the laser beam dwell time (diseussed in Chapter IV, Seetion D.2), is the 
eonduetive power loss from Equation (67) and is the radiative power loss from 

Equation (68). For the irradianee plots, the melt threshold will be identified by a solid, 
horizontal red line above whieh the laser will likely infliet signifieant damage to the 
target. 

The following plots present either irradianee or power-in-the-bueket data on the 
vertieal axis for a wavelength range of 0.5-5.0 pm, with the varying eolors representing 
the different values of the parameter being ehanged. Unless otherwise stated, the laser 
and atmosphere parameters eonform to the values listed in Tables 7 and 8. 
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1. Turbulence Study with a 5 km Target Range 

Figure 19 displays plots of the laser irradiance for the standard model, where the 

2 16 2/3 13 

Q parameter is varied from 10' m' (weak turbulenee, indicated by dark blue) to 10' 
m' (strong turbulence, indicated by light blue). As evidenced by the plot, turbulence 
corresponding to > 10”''^m“^''^ may prevent a successful engagement at a 5 km range, 

since the irradiance generally falls near, or below, the required melt threshold to inflict 
significant damage to the target. Furthermore, turbulence tends to affect the shorter 
wavelengths more than longer wavelengths (cf Equation 4), a trend reproduced by 
ANCHOR and visible in this figure. 


Standard Model Turbulence Study 

meitlhreshold _ 1e-13m''-2/3 1e-14m''-2/3 1e-15m^-2/3 16-16 m^>2/3 



1 2 3 4 5 

Wavelength (pm) 


Figure 19. Standard model turbulence study: Irradiance results indicate 
successful engagement at turbulences weaker than = 10~*‘'m“^^^. 

Figure 20 displays plots of the power-in-the-bucket for a bucket radius of 5 cm for 
the same study as Figure 19. Notice that moderate turbulence (C^ =10”'^m“^^^)does not 

decrease power-in-the-bucket, while stronger turbulence tends to increase overall on- 
target laser spot size. The moderate turbulence here does not increase the spot size 
enough to where it falls outside the bucket, and thus all of the laser power remains 
confined within it. Once the turbulence reaches extreme levels, however, the on-target 
spot size becomes larger than the bucket size, and thus the power-in-the-bucket drops. 
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standard Model Turbulence Study - PIB 

I e-13 m''-2/3 — e-14iii"-2/3 e-15ni*-2/3 1e-16 iti*-2/3 I 

10® I—-■---■---■---^-'-■-■-'-'-■-'— 



Wavelength (gm) 

Figure 20. Standard model turbulence study: Power-in-the-bucket results 
indicate significant loss of bucket power in the presence of strong 

turbulence. 

2. Variation of the Laser Beam Diameter at the Source with a 5 km 
Target Range 

The following study examines the trends in on-target laser irradiance and power- 
in-the-bucket associated with varying the laser beam diameter at the source at a target 
range of 5 km and in the presence of weak and strong turbulence. 

a. Weak Turbulence 

Figure 21 displays plots of the laser irradiance for the standard model, where the 
beam diameter at the source, and thus the beam director size D , has been varied from 
0.3 m (light blue) to 0.7 m (dark blue) within a weak turbulence (C^ =10“''^m”^^^) 

environment. In weak turbulence, increasing the beam director size shows marked 
improvement in the target irradiance. Furthermore, shorter wavelengths benefit less from 
the increase in the beam director size since turbulence and jitter contributions tend to 
dominate the time-averaged spot-size. 


62 














standard Model Source Laser Beam Diameter Study - Weak Turbulence 

melt threshold 0.3 m am 0.5 m 0.7 m 

10 ^ -—~:- 



1 2 3 4 5 


Wavelength (gm) 

Figure 21. Standard model weak turbulence beam director size study: Irradiance 
results indicate marked improvement in laser performance by 
increasing beam director size in a weakly turbulent environment. 


Figure 22 represents the power-in-the-bucket results for the same study as Figure 
21. Increasing the beam director size has little effect on the power-in-the-bucket in the 
presence of weak turbulence for shorter wavelengths; for longer wavelengths (>2 pm), 
the reduction of the spot size due to the larger beam director allows more of the power to 
be contained within the bucket. 


standard Model Source Laser Beam Diameter Study - Weak Turbulence - PIB 



Figure 22. Standard model weak turbulence beam director size study: Power-in- 
the-bucket results indicate little effect of beam director size on power 
delivered to bucket in weakly turbulent environment. 
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b. Strong Turbulence 

Figure 23 displays plots of the laser irradianee for the standard model, in 
which the laser beam diameter at the source, and thus the beam director size D , has been 
varied from 0.3 m to 0.7 m within a strong turbulence (C^ =10“'‘*m“^^^) environment. It 

appears that increasing the beam director size in this case has little effect on the 
laser performance, since the strong turbulence contribution dominates laser irradianee 
regardless of the size of the laser beam at the source. 


standard Model Source Laser Beam Diameter Study - Strong Turbulence 



Figure 23. Standard model strong turbulence beam director size study: 

Irradianee results indicate little effect of beam director size on laser 
performance in turbulent environment. 


Figure 24 displays plots of the power-in-the-bucket results for the same study as 
Figure 23. Increasing the beam director size has once again little effect on the power 
delivered within the bucket on the target within a strong turbulence environment. 
However, in comparison to Figure 22, the bucket power has markedly dropped for all 
wavelength windows; this occurs because in this more turbulent environment the spot 
size at the target is now larger than the bucket size regardless of the wavelength. 
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standard Model Source Laser Beam Diameter Study - Strong Turbulence - PIB 



Figure 24. Standard model strong turbulence beam director size study: Power- 
in-the-bucket results indicate little effect from beam director size on 
power delivered to a bucket in a turbulent environment. 

3. Meteorological Range (Haze) Study with a 5 km Target Range 

Figure 25 displays plots of the laser irradiance for the standard model, where the 
haze (meteorological range) has been varied from 5 km visibility (hazy day, indicated by 
light blue) to 23 km (clear day, indicated by dark blue). Strong haze (visibility < 5 km) 
prevents a successful engagement at a 5 km range, as the melt threshold is not likely to be 
met at such limited visibility. Furthermore, haze tends to affect shorter wavelengths more 
than longer wavelengths. As such, wavelengths greater than 3 pm experience less of a 
reduction in on-target irradiance than the 0.5 pm to 2.5 pm wavelength range. This is a 
known trend reproduced by these results [5]. 


standard Model Haze (Meteorological Range) Study 

^ - mell Ihreshold _ 5 km 10 km — 15 km 23 km 



1 2 3 4 5 

Wavelength (pm) 


Figure 25. Standard model meteorological range study: Irradiance results 
indicate successful target engagement at visibilities greater than 5 km. 
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Figure 26 displays plots of the power-in-the-bueket data for the same study as 
Figure 25. From the data, it appears that the reduction in power-in-the-bucket is 
proportional to the reduction in visibility. Instead of affecting the laser spot size, like 
turbulence does, haze instead attenuates the laser beam as it propagates to the target via 
atmospheric absorption and scattering. 



Figure 26. Standard model meteorological range study: Power-in-the-bucket 
results indicate overall reduction in power due to atmospheric 

attenuation. 


4. Precipitation Study with a 5 km Target Range 

Figure 27 displays plots of the laser irradiance for the standard model, where 
various amounts of precipitation have been introduced. Even the presence of light 
precipitation of 5.0 mm/hr rain prevents the likelihood of a successful engagement at a 5 
km range due to the attenuation of the laser power by the rain droplets. Irradiance 
estimates in heavier precipitation rates, while successfully calculated, fall below the scale 
of this plot at the 5 km target range. However, engagements at shorter ranges may still be 
viable, as will be seen later. 
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standard Model Rain Study - S km Target Range 


10 ” 


10 ° 


I 

I. 10^ 
0) 



10 = 


10 ’ 


melt threshold 
5.0 mm/hr Light Rain 


25.0 mm/hr Heavy Rain 
No Rain 



Wavelength (pm) 



Figure 27. Standard model precipitation study: Irradiance results indicate that 
even light rain prevents likelihood of a successful engagement at a 5 km 

target range. 


Figure 28 displays plots of the power-in-the-bucket data for the same study as 
Figure 27. Again, the results for heavier rain rates fall below this scale shown in the plot. 
A successful engagement is unlikely for even the slightest amount of precipitation at a 5 
km target range. 


standard Model Rain Study - 5 km Target Range - PIB 

25.0 mnVhr Heavy Rain 12.5 mm/hr Moderate Rain 5.0 mm/hr Light Rain 

No Rain 

106 . . -- 



1 2 3 4 5 

Wavelength (pm) 


Figure 28. Standard model precipitation study: Power-in-the-bucket results 

indicate significant power reduction in the presence of light rain at 5 km 

target range. 
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5. Beam Quality Study with a 5 km Target Range 

Figure 29 displays plots of the laser irradiance profiles for the standard model. In 
this figure the beam quality factor has been varied from 1 (perfect beam quality, 
indicated by dark blue) to 10 (poor beam quality, indicated by light blue). Longer laser 
wavelengths experience a greater reduction in performance with degraded beam quality. 
The reason is, once again, that longer wavelength spot sizes are dominated by vacuum 
diffraction (cf. Equation 45), where the beam quality plays a significant role. 


standard Model Beam Quality Study 

- meK threshold M''2=10 M''2=5 M'‘2=3 M''2=1 



1 2 3 4 5 


Wavelength (gm) 


Figure 29. Standard model beam quality study: Irradiance results indicate a 
greater reduction in laser performance at longer wavelengths with 
decreased beam quality. 


Figure 30 displays plots of the power-in-the-bucket data for the same study as 
Figure 29. Similar to the irradiance results, degradation in the beam quality causes a 
reduction in power-in-the-bucket for longer laser wavelengths, again due to diffraction. 
The increase in spot size reaches a point to where some of the laser power falls outside 
the bucket, correlating to the decreased power delivered to the desired target area. 
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standard Model Beam Quality Study - PIB 



Figure 30. Standard model beam quality study: Power-in-the-bucket results 
indieate greater reduction in bucket power at longer laser wavelengths 
for decreased beam quality. 

6. Platform Jitter Study with a 5 km Target Range 

Figure 31 displays plots of the laser irradiance profiles for the standard model, in 
which the platform jitter 0^^^^ has been varied from 0 prad (no jitter, indicated by dark 
blue) to 10 prad (extreme jitter, indicated by light blue). An increase in platform jitter 
results in a reduction of laser performance for all wavelengths; however, shorter 
wavelengths tend to be affected more. While the absolute contribution of jitter to the spot 
size is wavelength independent (cf. Equation 50), the relative contribution is greater for 
shorter wavelengths (cf. Equation 44). 


standard Model Platform Jitter Study 



Figure 31. Standard model platform jitter study: Irradiance results indicate an 
equal reduction in laser performance with increased platform jitter for 
all wavelengths, as expected. 
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Figure 32 displays plots of the power-in-the-bueket data for the same study as 
Figure 31. Similar to Figure 31, an increase in platform jitter increases the on-target laser 
spot size, which causes more laser power to fall outside the desired target area (the 
bucket), and thus less power is delivered within the bucket. 



Figure 32. Standard model platform jitter study: Power-in-the-bucket results 
indicate a reduction in laser performance with increased platform jitter 

for all wavelengths. 

7. Output Power Study with a 5 km Target Range 

Figure 33 displays plots of the laser irradiance results for the standard model, in 
which the laser output power P has been varied from 30 kW to 200 kW. From the 
results, it appears that increasing laser output power from 30 kW to 100 kW results in a 
marked performance increase of the laser. However, doubling the output power from 100 
kW to 200 kW has little benefit, as the increased laser power also increases thermal 
blooming, which causes a further reduction of on-axis intensity. Again, this is a known 
trend that is reproduced by the ANCHOR code [5]. 
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standard Model Output Power Study 
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Figure 33. Standard model output power study: Irradianee results indieate a 
marginal inerease in laser performanee when doubling laser output 
power from 100 kW to 200 kW. 


Figure 34 displays plots of the power-in-the-bucket data for the same study as 
Figure 33. While doubling the laser output power from 100 kW to 200 kW causes the 
effects of thermal blooming to be more pronounced (lowering the on-target laser 
intensity), power-in-the-bucket appears to be less affected by this phenomenon. While 
thermal blooming causes the on-target laser spot size to deform much like the depiction 
in Figure 6, the spot appears to remain generally confined within the desired on-target 
bucket, effectively delivering more power to the target at higher laser output powers and 
especially for shorter wavelengths. 



Figure 34. Standard model output power study: Power-in-the-bucket results 
indicate reasonable increase in laser performance by doubling output 
power from 100 kW to 200 kW. 
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c. 


STANDARD MODEL STUDIES AT VARIABLE TARGET RANGES 


The previous studies eoneerned the effeets of various atmospherie eonditions and 
laser properties on overall performanee while engaging a target at a fixed 5 km range. 
However, this approaeh is not wholly realistie as an adversary will rarely remain 
stationary in an operational environment. The following studies attempt to develop 
understanding of how atmospherie eonditions affeet laser lethality as the target eloses in 
on the laser souree. 

1. Target Range Variation Study 

Figure 35 displays plots of the laser irradianee results for the standard model, in 
whieh the target range has been varied from 1000 m (indieated by dark blue) to 7000 m 
(indieated by light blue). It appears that on a elear day with moderate turbulenee and a 10 
m/s erosswind, a target is likely to be sueeessfully engaged at ranges eloser than 7000 m 
for systems operating at the nominal ~1 pm SSL wavelength. 



Figure 35. Standard model target range study: Irradianee results indicate 
likelihood of a successful engagement at ranges closer than 7 km in 
clear weather conditions. 


Figure 36 displays plots of the power-in-the-bucket data for the same study as 
Figure 35. Similar to the conclusions drawn from Figure 35, as a target approaches the 
laser source, the laser lethality is increased. However, as the target closes in from 7000 m 
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to 1000 m, the power-in-the-bucket increases since less of the beam is attenuated by the 
atmosphere. At a range of -1000 m, nearly all of the power from the source falls within 
the bucket, since the spot size is relatively small even for the longer wavelengths. 



Figure 36. Standard model target range study: Power-in-the-bucket results 
indicate increase in laser performance as target closes in on laser 

source. 


2. Precipitation Target Range Studies 

The following irradiance and power-in-the-bucket studies evaluate laser 
performance at various target ranges and in the presence of precipitation of increasing 
severity. 

a. Drizzle 

Figure 37 displays plots of the laser irradiance results for the standard model in 
the presence of drizzle (2.5 mm/hr rain), with a variation in target proximity from 1000 m 
to 5000 m. While laser lethality in the presence of drizzle is marginal at best for a target 
range of 5 km, the on-target irradiance jumps by more than an order of magnitude as the 
target closes in to 3000 m, effectively making the weapon highly lethal against typical 
targets at this range. 
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standard Model - Drizzle (2.5 mm/hr) Target Range Study 
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Figure 37. Standard model drizzle target range study: Irradiance results indicate 
likelihood for successful engagement at ranges closer than 5 km. 


Figure 38 displays plots of power-in-the-bucket results for the same study as 
Figure 37. Similar to the results from Figure 37, as a target approaches the laser source, 
the power delivered to the bucket on the adversary increases for ah laser wavelengths. 
This increase in power-in-the-bucket, however, stems from the fact that the atmospheric 
precipitation attenuates laser power very effectively. The target has to be in closer 
proximity to the source to be hit by a substantial amount of laser power. 


standard Model - Drizzle Rain (2.5 mm/hr) Study - PI6 
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Figure 38. Standard model drizzle target range study: Power-in-the-bucket 
results indicate marked improvement in laser performance at target 
ranges closer than 5 km. 
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b. Light Rain 

Figure 39 displays plots of the laser irradiance results for the standard model in 
the presenee of light rain (5.0 mm/hr rain), with a variation in target proximity from 1000 
m to 5000 m. Similar to the results for 2.5 mm/hr rain, stronger precipitation reduces the 
on-target laser intensity by attenuating the laser power. The target needs to be within 
3000 m for the laser to effectively engage it. 



Figure 39. Standard model light rain target range study: Irradiance results 
indicate likelihood of successful target engagement at ranges closer 

than 3 km. 


Figure 40 displays plots of the power-in-the-bucket data for the same study as 
Figure 39. Similar to the power-in-the-bucket results for 2.5 mm/hr drizzle, a target range 
of less than 2000 m to the source results in an increase of power-in-the-bucket by almost 
an order of magnitude. 
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standard Model - Light Rain (5.0 mm/hr) Study - PiB 



Wavelength (pm) 


Figure 40. Standard model light rain target range study: Power-in-the-bucket 
results indicate marked increase in laser performance for target ranges 

closer than 3 km. 


c. Moderate Precipitation 

Figure 41 displays plots of the laser irradiance results for the standard model in 
the presence of moderate precipitation (12.5 mm/hr rain), in which the range to target has 
been varied from 1000 m to 4000 m. For shorter wavelengths the laser irradiance exceeds 
the melt threshold at a target range of 2 km. At a 1 km range, the on-target laser intensity 
exceeds the melt threshold by a large factor. 


standard Model - Moderate Rain (12.5 mm/hr) Target Range Study 
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Figure 41. Standard model moderate rain target range study: Irradiance results 
indicate likelihood of a successful engagement within a target range of 

2 km. 
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Figure 42 displays plots of the power-in-the-bueket data for the same study as Figure 
41. Similar to previous preeipitation results, the power-in-the-bucket inereases 
proportionally, and relatively equally for all wavelengths, to the target proximity from the 
source. However, the laser only truly becomes effective for target ranges closer than 2000 m. 


Standard Model - Moderate Rain (12.5 mm/hr) Target Range Study - PIB 



Wavelength (pm) 


Figure 42. Standard model moderate rain target range study: Power-in-the- 
bucket results indicate marked increase in laser performance for target 

ranges closer than 2 km. 


d. Heavy Precipitation 

Figure 43 displays plots of the laser irradiance results for the standard model in 
the presence of heavy precipitation (25.0 mm/hr rain), and as the target range is varied 
from 1000 m to 3000 m. The on-target laser irradiance exceeds the melt threshold at a 1 
km target range. Thus, even in heavy rain, a 100 kW-class HEL is a viable weapon 
against targets within 1000 m of the source. 


77 








Figure 43. Standard model heavy rain target range study: Irradianee results indieate 
likelihood of a sueeessful engagement at target ranges eloser than 1 km. 


Figure 44 displays plots of the power-in-the-bueket data for the same study as 
Figure 43. Similar to previous preeipitation results, power-in-the-bucket inereases 
proportionally with target proximity to the source. However, a 100 kW HEL only 
becomes effective as an offensive and defensive weapon in such heavy rains for target 
proximity of less than 1000 m. Otherwise, the actual power delivered to the target is 
much too low to inflict significant damage, although soft kills or blinding may still be 
achievable at longer ranges. 


Standard Model • Heavy Rain (25.0 mm/hr) Target Range Study • PIB 
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Figure 44. Standard model heavy rain target range study: Power-in-the-bucket 
results increase in laser performance for target ranges closer than 1 km. 
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3. Strong Haze Target Range Study 

Figure 45 displays plots of the laser irradiance results for the standard model in 
the presenee of strong haze (5 km meteorological range), within which the proximity to 
the target has been varied from 1000 m (indicated by dark blue) to 5000 m (indicated by 
light blue). In the presence of strong haze a target engagement is likely to be successful at 
a range of 5 km or less, as the melt threshold is exceeded at those ranges. The increase in 
on-target laser irradiance is dramatic; as the target approaches to within 3000 m, the laser 
intensity increases by more than an order of magnitude. 


Standard Model - Strong Haze (5 km VIS) Target Range Study 
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Figure 45. Standard model strong haze target range study: Irradiance results 
indicate successful engagement at target ranges closer than 5 km. 


Figure 46 displays plots of the power-in-the-bucket data for the same study as 
Figure 45. Closing the target in to the laser source results in more power delivered to the 
bucket. As the beam propagates a shorter distance, the hazy atmosphere attenuates the 
laser power less, and thus the on-target laser spot becomes more lethal. 
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Figure 46. Standard model strong haze target range study: Power-in-the-bueket 
results indicate increase in power delivered to the target for ranges 

closer than 5 km. 

4. Strong Turbulence Target Range Study 

Figure 47 displays plots of the laser irradiance results for the standard model in 
the presence of strong turbulence (C,^ =10”’‘*m~^^^), within which the proximity to target 

has been varied from 1000 m (indicated by dark blue) to 5000 m (indicated by light blue). 
A successful target engagement is likely for target ranges closer than 5000 m, as the melt 
threshold is exceeded at those ranges. Since the laser beam has to propagate a shorter 
distance to the target, the beam is less affected by the turbulent environment and the on- 
target spot size thus becomes smaller, while the laser irradiance becomes greater. 



Figure 47. Standard model strong turbulence target range study: Irradiance results 

indicate likelihood for successful engagement at target ranges closer than 5 km. 
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Figure 48 displays plots of the power-in-the-bueket data for the same study as 
Figure 47. It appears that as the target eloses in on the position of the laser souree, the 
power-in-the-bucket increases with the target proximity though less so at ranges closer 
than 2000 m. For these closer ranges, the effects of turbulence are so small that the 
power-in-the-bucket cannot increase with all of the available laser power already 
captured within it. 



Figure 48. Standard model strong turbulence target range study: Power-in-the- 
bucket results indicate increase in power delivered to target for ranges 

closer than 5 km. 


D. DWELL TIME STUDIES 

In addition to laser irradiance and power-in-the-bucket studies, the dwell time 
performance metric (discussed in Chapter IV, Section D.2) was also incorporated within 
the standard model to give additional insight on laser performance for various weather 
conditions and target ranges. 

I. Example Dwell Time Calculation 

Given a target with parameters as described in Section B, the required amount of 
time the laser spot needs to remain on target in order to inflict significant damage can be 
calculated by following the general procedure as discussed in Chapter IV, Section D.2. 
From Equation (67), the conductive power loss is calculated as 


81 








( 78 ) 


L 

and the radiative power loss from Equation (68) is given as 

'’«=a-£A(r:-r;) = 17W. (79) 

The total power loss to the surroundings is thus 

Pc+f’fi =2.4 kW. (80) 

In this example, the radiative power loss to the surroundings is rather negligible. 

Assuming that the laser delivers 50 kW of power to the target, the available 
melting power is ealeulated as =(50kW-2.4kW)»48 kW. Further assuming that 
the target is suffieiently damaged onee 0.1 L of its skin has been melted, the energy 
required to melt this volume «153 kJ. Finally, the dwell time required to melt the 

0.1 liters of aireraft skin ean be estimated by Equation (72) as 

^D=-^~3seo. (81) 

melt 

Of eourse, the aetual dwell times depend upon the speeifies of the target. The 
following studies aim to determine the dwell times required for a laser with the attributes 
listed in Table 7, exeept with a fixed 1.064 pm wavelength, to effeetively melt 0.1 liters 
of one millimeter-thiek aluminum. The speeified volume implies the assumption that the 
beam spot will wander while on-target in order to melt the requisite volume of material. 
The dwell times are determined for various haze, preeipitation and turbulenee eonditions 
and at different target ranges in order to fully quantify laser lethality. 

2. Dwell Time Study for Variable Visibility Conditions 

Figure 49 displays plots of the estimated dwell times required for this 1.064 pm 
wavelength laser to infliet damage to a target in various meteorologieal range eonditions. 
As evideneed by the plot, on elear days with visibility of 23 km, the required engagement 
dwell times fall below five seeonds with target ranges of 5 km or less. On hazy days with 
visibility of five km, the dwell times fall below five seeonds for target ranges of one 
thousand meters or eloser. Note that the 5 km visibility eurve diverges strongly at a range 


82 



of about three kilometers while the eurves eorresponding to less hazy eonditions do so at 
longer ranges. This is the result of the faet that at longer ranges, the eonduetive and, to a 
lesser extent the radiative, losses prevent the aluminum from reaehing its melting point 
regardless of the dwell time. 


Laser Dwell Time Study - Variable Haze (Meteorological Range) Conditions 
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Target Range (km) 


Figure 49. Laser dwell times for varying meteorologieal range eonditions; 

Results indieate reasonable engagement times for target ranges less than 
5 km (dependent on visibility). 


3. Dwell Time Study for Variable Precipitation Conditions 

Figure 50 displays plots of the estimated dwell times required for a laser to infliet 
damage to a target in various preeipitation eonditions. Even small amounts of 
preeipitation greatly affeet dwell times; however, target engagements within a few 
kilometers range remain viable. For heavy preeipitation, engagements are only viable for 
a target range of less than one kilometer. 
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Laser Dwell Time Study - Variable Precipitation Conditions 



Figure 50. Laser dwell times for varying precipitation conditions: Results 
indicate viable engagement times in rainy conditions at target ranges 

less than 1 km. 

4. Dwell Time Study for Variable Turbulence Conditions 

Figure 51 displays plots of the estimated dwell times required for a laser to inflict 
damage to a target in various turbulence conditions. Dwell times remain reasonable even 
in a very turbulent atmosphere for target ranges of five kilometers or less. 


Laser Dwell Time Study - Variable Turbulence Conditions 
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Figure 51. Laser dwell times for varying turbulence conditions: Results indicate 
reasonable target engagement times for ranges less than 5 km, even in 
the presence of turbulence. 
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VII. CONCLUSIONS 


A. CONCLUSIONS FOR 5 KM TARGET RANGE 

According to the results for the standard model studies at a 5 km target range, the 
100 kW-class HEL is most likely to successfully engage slow-moving targets (such as 
light aircraft, UAVs, etc.) only in mild atmospheric conditions. Specifically, this 
requirement is met if atmospherie turbulenee is weak to moderate (C^ =10“''’m“^^^ to 

lO^^^m^^^^), the meteorological range is greater than 10 km visibility, and no appreciable 
preeipitation is present. 

Furthermore, it appears that large beam director sizes increase on-target laser 
irradianee in only weak turbulenee eonditions. It is quite possible that moderate to strong 
turbulence would be present on some days while in a maritime environment. So, laser 
weapon systems installed on naval vessels or aircraft could be designed in a more 
compact manner with smaller beam directors, since larger beam directors provide little 
added benefit in sueh conditions. 

B. CONCLUSIONS FOR TARGET RANGES CLOSER THAN 5 KM 

According to the results for the standard model studies at variable target ranges, 
the 100 kW-class HEL can be effectively utilized in an offensive and defensive capability 
against slow-moving targets (such as light aircraft, UAVs, etc.) in a greater variety of 
atmospheric conditions as the target range decreases below five kilometers. Since the 
various eontributions to on-target laser irradianee (diffraetion, jitter, turbulence, 
extinetion and thermal blooming) are strongly dependent on target range, these 
eontributions eonspire to greatly inerease laser irradianee at shorter distanees. 

Furthermore, these effects also aid in the total power delivered to the desired 
target area, the bucket. The closer a target is to the laser source, the lesser impact that 
atmospheric attenuation and beam dispersion have on the laser beam as it propagates to 
the objective. From these results, 100 kW-class HELs could be effeetively employed 
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against targets in hazy, turbulent, and rainy conditions as long as the target is no more 
than three kilometers from the laser source, depending on the severity of the weather. 


C. CONCLUSIONS FOR VARYING DWELL TIME STUDIES 

According to the dwell time studies in clear weather, the 100 kW-class HEL can 
successfully engage and damage a slow-moving target from a range of five kilometers or 
closer within a five-second dwell time. For target ranges from 5 km to 10 km, longer 
laser dwell times (~15 seconds) are required to inflict a similar amount of damage. For 
target engagements at longer ranges and in poor weather conditions, required laser dwell 
times rise to greater than 10 seconds, and in some cases become unattainable (in the case 
of heavy rain, for example). However, for hazy, turbulent and rainy days, laser dwell 
times are still reasonable at shorter ranges (< 3 km). 
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